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ABSTRACT 


A unified theory of alignment is developed which is 
founded on the two considerations common to all alignment, pro- 
cedures: coordinate frames and measurements. All alignment 
procedures are viewed as a class of non- dimensional , vector - 
direction-indication problems which effect a transfer of orien- 
tation between two coordinate frames. Part of the present the- 
sis is that the differences among Various approaches to the 
solution of the moving base alignment problem are superficial. 

A discussion of the essential requirements for alignment is 
presented in the unified format. The roles of instrument errors 
and base motion as fundamental limitations on the alignment ac- 
curacy are developed in detail. All performance indices are 
chosen as angular quantities which are invariant under coordinate 
transformation. This facilitates a comparison and evaluation of 
techniques independent of the details of system mechanization. 

The necessary conversion of magnitude -sensitive instrument 
uncertainties to their angular equivalents is presented in detail. 

Purely geometric considerations yield a closed form 
solution for the sensitivity of platform misalignment to the angle 
between measured vectors. It is also shown that the magnitude 
of a rotation specified, by the arguments of the principal direction 
cosines is independent of the skewness of the rotated coordinate 
frame. The magnitude of a rotation specified by Euler angles is 
found to be almost completely independent of the order of the ro- 
tations. The very important distinctions among "measuring a 
vector, " "measuring a vector's magnitude, " and "measuring a 
vector's direction" are quantified as a function of instrument 
error parameters and the orientation of the vector with respect 
to the instruments. 

Application of the unified analysis is made to the following 
specific alignment techniques: Vertical Indication; Gyrocom- 
passing; Star Tracking; Fix Monitored Azimuth; Vector Matching; 
Gimbal Angle Matching; and Optical Comparison. It is shown 
for a multiple -system configuration that simultaneous measure- 
ments of a single vector by multiple systems do not constitute 
multiple measurements of a common alignment parameter. An 
additional measurement is introduced which takes advantage of the 
multiple system configuration and improves the alignment of 
each system. The extent to which operational considerations act 
as a pre-filter on the selection of an alignment scheme is dis- 
cussed relative to the carrier-borne alignment of aircraft navi- 
gation systems. 
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ALIGNMENT OF INERTIAL SYSTEMS 
ON A MOVING BASE * 

By Arthur H. Lipton 
Electronics Research Center 
Cambridge, Massachusetts 


CHAPTER I 
INTRODUCTION 


The preparation of an inertial guidance system for 
any mission requires: 

1. Establishment of a measurement reference 
frame through initial alignment of the inertial 
measuring unit (IMU), and 

2. Establishment of initial conditions on 
velocity and position of the system in this frame. 

The initialization of a system while its carrying vehicle is 
stationary relative to the earth has been the subject of a 
continuing effort to achieve greater accuracies in less time. 

This effort has met with considerable success in the past, and, 

^ stisfc 

in fact, is currently showing significant progress. 

There has also been a recurring requirement to per- 
form this initialization when the carrying vehicle is non- stationary 
relative to the earth. The concept of initialization- in- motion 
is particularly attractive in weapon system technology. Several 
examples of vehicles which may be used to transport an inertial 
system together with representative systems which specify 
(or specified) initialization while in motion are an aircraft 
(Skybolt, Hound Dog), a train (Minuteman), an aircraft carrier 
(F-111B) or a submarine (Polaris). 

* Submitted to the Massachusetts Institute of Technology, Department of 
Aeronautics and Astronautics in August 1966 in partial fulfillment of 
the requirements for the degree of Doctor of Science in Instrumentation. 

** Superscript numerals refer to numbered references in the 
bibliography. 



Unfortunately, success in initialization on a moving 
base has been much more limited than that on a stationary base. 
The initial alignment portion of the complete initialization task 
represents a particularly hard core problem. It is to this 
initial alignment problem that this thesis is addressed, with 
special emphasis on the problems of aligning on a moving base. 

Alignment, for our purposes, is defined as deter- 
mining the angular orientation of a set of fiducial axes fixed in 
an IMU with respect to an arbitrarily chosen set of reference 
axes. This determination is either an analytic evaluation of the 
relative angular orientation between two coordinate frames or a 
procedure which seeks to achieve angular- coincidence between two 
coordinate frames. The reference axes to which the IMU is 
aligned frequently represent a traditional geographic or geocen- 
tric coordinate frame. In several moving base alignment tech- 
niques, however, the reference axes may well be the fiducial axes 
of a second IMU. Regardless of the particular frames involved, 
the error quantity (misalignment) represents an undesired angular 
discrepancy between the IMU' s fiducial axes and the reference 
coordinate axes. The most general description of misalignment 
is the "whole angle rotation" described in Chapter II. A familiar 
and important specialization of this description is the vertical 
and azimuth specification of an IMU' s misalignment relative to 
a geographic reference frame. 

Most alignments depend, to some degree, on measure- 

2 

ments of the IMU' s inertial sensors : accelerometers and gyros. 
Since these sensors are required for guidance and navigation 
subsequent to initialization, it is desirable to utilize them for 
alignment. The accelerometers and gyros are able to effect 
alignment relative to the earth by identifying the earth's gravity 



(g) and daily rotation rate (u j-g) vectors. This identification is 
straightforward from a base fixed relative to the earth. Motion 
of a vehicle with respect to the earth, however, effectively de- 
couples the system from these significant inputs' by adding iner- 
tially indistinguishable vehicle accelerations to g and angular 
vehicle rates to uJ When this decoupling due to base-motion 
becomes severe enough to prohibit alignment by standard fixed- 
base techniques, we have a moving base alignment problem. As 
used in this thesis, then, a moving base connotes a vehicle which 
undergoes sufficient translation and rotation relative to the earth 
to seriously affect the alignment problem. 

An examination of the literature ^ reveals two 
basic engineering approaches to the solution of the moving base 
alignment problem. The first is an extension and/or modification 
of the techniques used in a stationary situation. The extensions 
and/or modifications effectively serve to recouple the system to 
the earth by processing (i. e. filtering, compensating) the moving- 
base inertial measurements in a manner which obtains an equiv- 
alent fixed-base inertial measurement. Often as not, this requires 
additional non- inertial measurements with respect to the earth. 

The second approach is the development of a group of 
special techniques. Some of these techniques still rely upon 
measurements of the systems' inertial sensors, but the measured 
quantities no longer have a unique a priori meaning relative to the 
earth. In this case, some means of calibrating the measured 
quantity must be provided, usually in the form of a pre- aligned 
inertial measurement unit carried on board the same vehicle. 

Other special techniques introduce non- inertial alignment 
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measurements as a means of circumventing the base motion 
disturbance of inertial quantities. Some of these non- inertial 
measurements identify quantities which have an a priori sig- 
nificance relative to an inertial reference frame (e.g. , a star 
line- of- sight) or even to the earth (e.g., a range measurement 
relative to an earth- fixed point). Still other non- inertial 
measurements are merely a means of transferring the alignment 
of a pre- aligned on-board IMU to the inertial system of an 
aircraft or missile. 

This thesis emphasizes the homogeneity among 
the superficially distinct approaches to the moving base alignment 
problem. Among other things, this facilitates the evaluation 
and comparison of the more significant techniques. The funda- 
mental limitations of moving base alignment are viewed as 
1) the errors due to the measurement devices and 2) the errors 
due to base motion interference with the measurement paths. In 
this context, base-motion interference relates not only to gross 
motions of the vehicle relative to the earth but also to the 
non-rigid behavior of the vehicle structure. That these fundamen- 
tal limitations may be described in a uniform fashion regardless 
of the measurement technique is demonstrated in the thesis. The 
relative contributions of sensor error and base- motion to a 
composite alignment accuracy are developed, providing a basis 
for studying the trade-off between instrument errors and disturb- 
ing base- motions. 

Chapter II of this thesis is devoted to the development 
of a unified theory of alignment. This theory is founded on the 
two considerations common to all alignment procedures: 
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coordinate frames and measurements. All alignment procedures 
are viewed as the transfer of alignment between two coordinate 
frames by means of measurements on at least two non-collinear 
vectors. Individual techniques are seen to vary only in the 
mechanization of particular measurements. Coordinate frame 
rotations are treated as whole angle rotations, the latter being 
developed in terms of all the more usual rotational parameters. 
Measurements are discussed from the viewpoint of single- 
degree- of-freedom sensors attempting to locate the direction 
of a vector input. Both magnitude sensitive and angle sensitive 
instruments are considered. Base motion is treated as 
causing an angular deflection of the measured vector from a 
desired nominal direction. The magnitude of the total platform 
misalignment resulting from independent measurements of.two 
non-collinear vectors is calculated as a function of the angular 

I 

uncertainties in determining the nominal direction of each of 
the vectors. 

Chapter III applies the analysis of Chapter II to the 
following specific alignment techniques: Vertical Indication, 
Gyrocompassing, Star Tracking, Fix Monitored Azimuth, Vector 
Matching, Gimbal Angle Matching and Optical Comparison. The 
mathematical methodology not only provides a common basis for 
comparing techniques but also is seen in certain instances to be 
a simpler way of treating old problems (e. g. , Gimbal Angle 
Matching). Limitations on the use of each technique are discussed 
as an aid in the appropriate selection of measurement techniques 
for individual situations. Individual conclusions are necessarily 
limited in scope but taken together encompass the spectrum of 
moving base problems associated with acceptable carrying vehicles. 
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Chapter IV discusses two unique operational situations. 
The first is that of multiple missile systems carried aboard the 
same vehicle. The question addressed is whether the simultaneous 
alignment of all the systems is advantageous from the standpoint 
of improved individual system alignments. The second situation, 
that of the carrier-borne alignment of aircraft navigation systems, 
demonstrates how operational considerations can be a very strong 
pre-filter on the selection of a particular technique. 

A summary and recommendations are presented in 

Chapter V. 

Professor Walter Wrigley, the chairman of my thesis 
committee, has been a source of encouragement and helpful 
suggestions throughout my doctoral program. Detailed discussions 
with him and with Professors Y. T. Li, W.R. Markey and 
J. E. Potter have contributed greatly to this thesis and are sin- 
cerely appreciated. 

I thank Mr. Kenneth Fertig for the encouragement and 
training received while under his supervision and for his continued 
interest in the author's progress. I am also grateful to my col- 
leagues at the MIT Instrumentation .Laboratory and the NASA 
Electronics Research Center for their many helpful suggestions 
and discussions. I acknowledge with gratitude the understanding 
of my supervisors at NASA in providing conditions necessary for 
completion of this thesis. 
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CHAPTER II 


A UNIFIED THEORY OF ALIGNMENT 


2. 1 I ntroduction 

The literature devoted to the alignment of inertial systems on 
a moving base (for example,- References 1, 2, 3, 4, 13, 24, 25, 27, 

30, 32) is primarily concerned with the analysis of distinct approaches 
and techniques. In fact, seven specific alignment techniques are 
discussed in detail in Chapter III of this thesis. Yet no matter how 
different the various approaches to the solution of the moving base 
alignment problem may appear, this difference is superficial. There 
exists a large measure of homogeneity among the various approaches 
which is not evident from a survey of the literature. Properly em- 
phasized, this homogeneity facilitates a comparison and evaluation of 
techniques which a preoccupation with disparity might otherwise mask. 
Hence the development and presentation of this unified theory of 
alignment. 

When reduced to its elementary form, the alignment problem is 
seen to be concerned with three major considerations: l)the consi- 
deration of coordinate frames and coordinate frame transformations; 

2) the consideration of measurements made by sensors integral to a 
particular inertial system; and 3) the consideration of information 
transferred from an external sensor or system. 

The organization of this chapter is such that these considerations 
are first discussed separately. This is primarily to minimize the 
framework necessary for analysis of the problem. It is not to imply 
that the three (3) considerations are completely independent, for as 
measurements are made to establish the alignment of coordinate 
frames, so does the alignment of coordinate frames affect the ability 
of sensors to make accurate measurements. Dependencies such as 
this are developed as part of the uniform error analysis in the latter 
portion of the chapter. 
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2.2 Coordina te Frame Cons iderations 

The definition of alignment requires only two (2) coordi- 
nate frames: a set of fiducial axes fixed in an inertial measuring 
unit (IMU) and an arbitrary set of reference axes. Yet the phy- 
sical realization of an alignment procedure frequently requires 
the consideration of many more than two (2) coordinate frames. 

By minimizing the required number, one achieves a significant 
simplification in the associated transformations and manipulations. 

It is possible to identify four reference coordinate frames 
which constitute the absolute minimum necessary to describe 
every moving base alignment procedure known to this author. 

This minimum results from considering the stable member of an 
IMU as simply a rigid body. In order to permit evaluation of 
the alignment (i. e. , orientation or attitude) of this stable mem- 
ber, one need only visualize implanting a reference triad in the 
stable member. Except for certain purposes of mathematical 
convenience, the location of this reference set does not have to 
have a special relationship to instrument axes, computational 
axes, or the navigational reference frame. It will, however, or- 
dinarily be chosen as an orthogonal set. 

Figure 2. 1 shows the functional relationships for the re- 
quired reference coordinate frames. In addition, because the 
alignment process is always based on a set of measurements 
and sometimes aided by a transfer of external information, the 
alternate measurement paths required to describe the complete 
alignment problem are also indicated. Of the two coordinate 
frames required for the definition of alignment, the fiducial axes 
of an IMU form the dependent reference frame of Figure 2. 1 and 
the arbitrary set of reference axes form the independent refer- 
ence frame. 

2.2.1 Independ e nt Reference Frame 

The independent reference frame is a purely ana- 
lytical coordinate frame. As such it is the only required coor- 
dinate frame which is not considered rigidly attached to the stable 
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Figure 2-1 Functional Diagram of Coordinate Frames and Measurement Paths Required 

to Describe the Moving Base Alignment Problem 



member of an IMU. It represents the ideal orientation of any of its 
physically mechanized counterparts and is therefore the ultimate 
gauge of how well a system has been aligned. The orientation chosen 
for the independent reference frame usually will be closely related to 
the geometry of the space chosen for navigation and guidance of the 
carrying vehicle. The traditional independent reference frames, such 

as geographic, geocentric, and inertial are described in detail in the 

1>t . 5, 6, 7, 8 

literature 

2.2.2 Dependent R eference Frajrie 

Each of the dependent reference frames indicated in 
Figure 2. 1 is associated with one IMU which is to be aligned. The 
multiplicity recognizes that a single vehicle might carry several sys- 
tems which require alignment. In the event that these systems are 
different, the dependent reference frames may be non-identical. Con- 
versely, identical dependent frames are expected for identical systems. 
Chapter III discusses several techniques whereby each dependent frame 
is aligned individually. Chapter IV considers the simultaneous align- 
ment of several dependent frames. 

The dependent frames are reference frames only inso- 
far as they serve as a measurement reference for the IMU during a 
mission. They are dependent because they must be aligned with re- 
spect to some other reference system. In particular, it is desirable 
to align the dependent system by a series of measurements made with 
respect to the independent reference frame. This requires informa- 
tion which can be sensed by the instrument associated with the depen- 
dent reference frame to have a unique geometrical relationship to the 
independent reference frame. For example, the earth's gravity vector 
(g) and angular velocity vector ( w are well defined in the geographic 
frame. If their orientation with respect to the dependent frame can be 
established by measurements performed in the dependent frame and 
compared to the known orientation in the geographic frame, the align- 
ment of the dependent frame to the independent frame follows directly. 
This particular procedure is known as gyrocompassing. It depends 
upon the capability to measure g and u ^ accurately in the dependent 
reference frame. 



2.2.3 Intermediate Reference Frame 


When measurements suitable for alignment can- 
not be made directly with respect to the independent reference 
frame, an indirect series of measurements may be used. For 
this purpose a "substitute" independent reference frame, i. e. 
one which "remembers" the independent reference frame, is 
provided. This is the intermediate reference frame of Figure 
2. 1. The use of active intermediate reference frames for mo- 
ving base alignment is perhaps the most significant dividing 
point between fixed base and moving base alignment procedures. 
Although intermediate reference frames are often necessary in 
both fixed and moving base alignment to circumvent physical 
disturbance of the measurement, the fixed base procedure will 
usually employ a static type of intermediate reference. In 
many proposed moving base alignment techniques, there is made 
available on board the vehicle an IMU which has previously been 
aligned with respect to an independent reference frame and to 
which it is desired to align one or more additional platforms. 

In other words, it is desired to transfer the alignment of the 
intermediate reference frame to the dependent frame. In this 
case the alignment capability is measured with respect to the 
intermediate reference frame while the total alignment accuracy 
must include the ability of the intermediate system to remember 
the orientation obtained from the independent reference frame. 

A transfer alignment of the type considered above is some- 
times called a Master-Slave alignment. This name derives from 
mechanizations where the orientation of the dependent system 

(slave) is physically slaved to match the orientation of the inter- 
3 

mediate system , but is now popularly applied to any transfer 
5 4 

method Gulland has discussed a naval application involving 
transfer of the SINS coordinate reference to a remote fire- 

3 

control device aboard the same vessel. Nauman and Oestreich 
5 

and O'Donnell indicate the interest in a transfer of coordinates 
from an airborne auto- navigator to a remote fire-control device 
aboard the same aircraft. 



2.2.4 Transfer Reference Frame 

When it is not possible to make suitable measure- 
ments between dependent and intermediate reference systems, a 
transfer reference frame (Figure 2. 1) may be used. The proposed 
carrier-borne alignment of the F-111B navigation equipment* ® 
employs this technique. The transfer reference frame is an active 
inertial platform which is first aligned to the intermediate system. 
Then it is physically transported to where its alignment can be 
measured by and transferred to the dependent reference frame. 
Although required to remember its orientation for a much shorter 
time period than the intermediate system, the transfer system and 
errors introduced by the alignment procedure are similar to those 
considered for the intermediate system discussed in Section 2. 2. 3. 

2. 3 Measurement Consideratio ns 

2. 3. 1 Vectorial Nature of the Measurements 

All alignment techniques depend on the ability to ac- 
curately measure naturally existing or artificially introduced 
quantities. The point of greatest significance to the unification 
of alignment studies is that regardless of the fact that the mea- 
surement of these quantities may proceed by inertial, optical, 
electromagnetic or other means, the various measurements 
have a common characteristic of singular importance: 

Every quantity measured for the purpose of alignment is 
vectorial in nature. 

The vectorial nature to which we refer is the particular concept 
of a vector as a geometric entity having characteristics of direc- 
tion and magnitude. Although every measured quantity relevant 
to alignment may, in fact, be properly described as a vector, the 
determination of just its direction supplies conditions which are 
mathematically sufficient to determine alignment. The determination 
of magnitude is not a mathematically necessary condition but it may 
be a practical necessity depending upon the nature of the measured 
quantity and the sensors employed. When magnitude measure- 
ments are made, they are for the purpose of determining the 
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directional parameters of the input quantity. This concept is 
sufficiently important to warrant amplification by example. 

Consider, first, the fixed base gyrocompassing 
mode of alignment. Here the system attempts to establish its 
orientation by determining the vectorial direction of each of 
two physical quantities: gravity (g) and the earth's daily ro- 
tation (u . ). Each of these quantities is a vector. The instru- 
mented measurements (using accelerometers and gyros respec- 
tively) determine the vector directions in one of two ways. 

Either they make the best possible estimate of magnitu de com- 
ponents of the quantity or they seek the null magni tude plane 
of the quantity. This presents a choice of four possible measure- 
ment combinations, but in each case the full vector nature of 
gravity and earthrate is being utilized. 

Optical measurements provide an example of an 
alignment measurement which is not a complete vector mecha- 
nization but whose basis is indeed vectorial. Optical systems 
(such as auto- collimators) are concerned with the line of sight 
from the instrument to an object. This line of sight may be 
prbperly considered to represent only the directional aspect of 
the position vector from the instrument to the object, where 
in order to rpeasure rotation of the object about some axis the 
position vector must terminate in a point of the body off the rota- 
tional axis. Angular measurements can determine the direction 
of this position vector (and hence the alignment between object 
and instrument) without recourse to magnitude measurements. 

In optics, magnitude considerations relate to the energy of the 
electromagnetic wave. In this sense, the magnitude of the po- 
sition vector does weigh heavily on the practical application of 
optical instruments. Other parameters of the wave, however, 
carry the intelligence related to direction finding. 

It may now be seen that examining the alignment 
problem in terms of a vector measurement problem permits a 



clear understanding of the limitations on the process of moving base 
alignment implied by the measurement capabilities. It is primarily 
in the choice of the vectors and the corresponding form of the 
measurement that the techniques discussed in the literature 
differ . 

2. 3. 2 Two Vector Requir ement 

The determination of the angular orientation of 
an inertial measuring unit by vector measurement requires at 
a minimum two non-collinear* vectors whose orientation is known 
with respect to the desired system attitude. This may be proved 
simply by both mathematical and physical arguments. 

Consider a misalignment between the dependent 
and independent frames of Figure 2-2. A relationship between 
the observations of the vector V in each frame is 

V 1 = V d 2. 3-1 

where the vector superscripts indicate coordinatization in the 
respective frames and the rotation matrix carries the d- 
frame into the i-frame. Visualize the misalignment as being 
developed on an incremental basis beginning with the two frames 
coincident and having V remain constant in the i-frame. Then 
we may write 

AV 1 = R* [AV dj -R d AR 1 . V d ] 

d l a 


from which 


= 0 

- AV d = R d AR 1 V d 
i d 


2. 3-2 


* In the general case, one merely requires a pair of skew vec- 
tors. Because a center of measurement is defined in the IMU, 
the measured vectors must intersect and define a common plane. 
Hence the requirement for non-collinear vectors. 
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Figure 2-2 Evaluating Misalignment (d0) Between Independent (i) 
and Dependent (d) Frames by Measurement of a Vector (V) 



Has 3 degrees of freedom about Has one degree of freedom about 

any three independent axes. measured vector V. 



Orientation fixed with respect to two vectors, 
0 degrees of freedom. 


Figure 2-3 A Physical Argument to Establish the Two-Vector Requirement 
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or 



R d AR 1 , V d 
1 d 


2. 3-3 


where we choose t- 0 as the instant when the d-frame and i- 
frame are coincident. We ask now whether these observa- 
tions of the single vector V are sufficient to evaluate uniquely 
the misalignment between the frames. 

Under suitable restrictions (which have been 
11 12 13 

discussed by many authors ' ' ) the misalignment may 

be represented by a vector angle d0 . The direction of d 9 gives 

the axis about which the misalignment was generated, while 

the magnitude of d# specifies the amount of rotation about this 

axis. Using the vector angle concept, an identification of the 

right-hand side of Equation 2. 3-3 in terms of a vector cross 
14 

product may be made. 


R d ARjj =>d6x 


2. 3-4 


Substituting 2. 3-4 into 2. 3-3 gives the desired expression, 
coordinatized in the d-frame 

V d j t!<) - V d = dflxV 2.3-5 

The non- simultaneous measurements of 2. 3-5 may be replaced 
by equivalent simultaneous observations of the same vector 
from different frames, since 





t=0 



t=0 


= V 1 2. 3-6 

— d 1 

The term V | q in Equation 2. 3-5 is, therefore, a represen- 
tation in the d-frame of V 1 . While the substitution of 2. 3-6 into 
2. 3-5 is rigorously proper, the resulting equation is subject to 
mis-interpretation and therefore is not written explicitly. 
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The implication of Equation 2. 3-5 is that only 
the component of d0 perpendicular to V can be evaluated by mea- 
surements of V. Any component of d0 which is parallel to V re- 
mains indeterminate after a measurement of V. The evaluation 
of this parallel component requires measurement of information 
which is linearly independent of V in three- space, hence has some 
component orthogonal to V. Because the additional information 
must here be described as a vector, the requirement for linear 
independence is met by selecting the second vector to be non- 
collinear with V. The component of the second vector which is 
normal to V contains the required intelligence; the component 
parallel to V supplies redundant information which may or may 
not be used to improve the estimate of the d0 component perpen- 
dicular to V. Note in this connection that a null- magnitude com- 
ponent of a vector is equally as signficant for direction indication 
as a non-zero projection. However, the greatest, sensitivity for 
evaluation of the parallel component of d0 obviously corresponds 
to choosing the second vector perpendicular to V. 

A simple physical argument establishing the 
two vector requirement proceeds as follows. The stable member 
of an IMU, considered as a rigid body, has three degrees of ro- 
tational freedom (Figure 2-3 a). The determination of the body's 
orientation with respect to one measured vector is entirely ana- 
logous to fixing that vector in the body, much as an axle (Figure 
2-3 b). The original three degrees of freedom are thus reduced 
by two, leaving only the single degree of rotational freedom about 
the originally measured vector unspecified. The required non- 
collinear vector (Figure 2-3 c), by supplying a measurement ortho- 
gonal to V, fixes the orientation of the rigid body. 

The visualization of a measured vector as an axis in the 
body, when specialized to the case of a coordinate set "straddling" 
a vector input, is of central importance in developments to follow. 



It is emphasized that by measuring two non-collinear 
vectors it is possible only to uniquely specify the orientation of the 
dependent reference frame with respect to the vectors. The orienta- 
tion of the vectors must be known with respect to the desired system 
attitude (as represented in the independent, intermediate or transfer 
reference frames) in order to complete the measurement portion of 
the alignment. In most cases it is interesting to make a detailed ex- 
amination of the accuracy with which a system assumes a new orien- 
tation from a measured orientation only when the dynamics of this 
process affect the measurement problem itself. 

An important violation of the non-collinear vector 
requirement occurs when a gyrocompassing system is operated in the 
polar regions. The often discussed singular behavior of such systems 
has a simple physical basis in the near-collinearity of g and w 
Precisely at the poles, where g and to become collinear, the vec- 
tors no longer provide sufficient information to determine direction 
uniquely. 

Since each measured vector provides two of the three 
pieces of information required for alignment, a mathematical over- 
specification of alignment results from the use of two non-collinear 
vectors. This leads quite naturally to the possibility of measuring 
more than two vectors in order to take advantage of this redundancy. 

In fact, measuring a total of three vectors, no two of which are colli- 
near, would provide exactly twice the conditions mathematically suf- 
ficient to determine alignment; however, this approach has little 
practical appear at present because of the difficulty involved in instru- 
menting even the two required vector measurements. The specific 
problem of employing more than two vectors is not pursued in this 
thesis, but much of the forthcoming analysis is applicable to this concept. 

In general, the requirement for two non-collinear 
vectors cannot be met by two vectors of the same type if the measure- 
ments are made simultaneously. The exception to this is represented 
by the direction vectors involved in optical measurements. This is 
because two physical vectors with the same units of measure will add 



vectorially and appear to a measuring instrument as a single re- 
sultant vector . For example, although u) and another non- 
collinear angular velocity vector may exist at the same time, 
they appear as just a single vector to the gyroscopic sensors. 

For a simultaneous measurement we must therefore require 
(with the exception noted above) that the two vectors employed be 
distinguishable on a dimensional basis in addition to being non- 
collinear . 

Non- simultaneous measurements make it pos- 
sible not only to use two different vectors of the same type 
(providing they do not exist simultaneously) but to have a single 
physical vector satisfy the two vector requirement. This requires, 
however, that the direction of the vector change as a function of 
time in both the reference frame to which we are aligning and 
the dependent frame being aligned. By way of explanation, con- 
sider a different viewpoint of fixed base gyrocompassing. If a 
system's gyros are used to maintain an inertially fixed orien- 
tation, observations of the local gravity vector display a uniform 
change of direction relative to the system. With the aid of a clock, 
these observations of a single vector provide sufficient informa- 
tion for alignment relative to an inertial independent reference 
frame. The important point is that the physical vector (g) has re- 
mained constant; only a carefully calibrated apparent direction 
change has occurred. It is this apparent change which permits 
alignment. 

Vector matching (Section 3. 6) provides an exam- 
ple where non- simultaneous measurements of similar vectors 
permit alignment. The master and slave systems are both aligned 
to "g with the vehicle in a non- maneuvering condition. The systems 
are then placed in a vertical-hold mode while the vehicle maneuvers 
in order to generate a significant horizontal acceleration. Although 
the total specific force vector now lies somewhere between g and 
the horizontal, our "calibration" of g at a prior time period allows 
us to complete the alignment based upon the horizontal acceleration 
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only. Had we not separated the measurements in time, the 
components existing at one point in time would be indistinguish- 
able. 

It is crucial in these mechanizations that a real 
or apparent change of vector direction occur in both the depen- 
dent and independent frames. If the vector orientation were to 
remain fixed with respect to either coordinate frame, there 
would be an unresolved ambiguity in that frame's orientation 
about the vector. An ambiguity associated with either frame, 
is, of course, an alignment uncertainty between the two frames. 
In order for the non- simultaneous vector technique to be con- 
sidered, the reference frame must be capable of monitoring 
the vector's direction change very accurately. It should be 
obvious now that the actual change of vector direction may be 
accomplished either by rotating the reference frame and de- 
pendent frame (called "slewing") with respect to a spatially 
fixed vector or by spatially fixing the coordinate frames and 
turning the vector. 

The foregoing discussion of the measurement 
problem has indicated the principal choices available for meeting 
the minimum vector measurement requirement but has not been 
concerned with the quantities the input vector represents. If 
one counts the combination of ways available to meet the two 
vector requirement while also distinguishing the vector quanti- 
ties on the basis of dimension (i. e. , position, velocity, accele- 
ration . . . ), the selection of an appropriate solution of the 
alignment problem is already seen to be a formidable task. 

This selection is further complicated by the external measuring 
equipment and information transfers depicted in Figure 2.1. 
Separate from the measurements made by sensors integral to 
the inertial system under consideration, information may also 
be transferred to the system from a remote sensor or system. 
This information is considered separately from primary mea- 
surements not only because its role is usually auxiliary but also 



because the separation between measurement points allows the 
information link to be corrupted by base motion interference effects. 
The important point to recognize* and which renders the general 
analysis of moving base alignment tractable, is that because every 
input is a vector quantity, the error analysis may be treated by 
simply relating fixed directional indicating capabilities to a resul- 
tant alignment accuracy. 

The inaccuracy of the measurement processes, has 
two major sources: the inherent errors of the measuring instru- 
ments and the fact that all measurement paths are subject to 
interference from base motion. The success (or failure) of 
present alignment procedures can be said to depend largely on 
how well base motion is either 1) circumvented, 2) measured, 
or 3) predicted and filtered, since the instrument quality is usually 
taken as a fixed item. The approach of this thesis is to show the 
tradeoff required between the two error sources in order to achieve 
a desired alignment accuracy. 

Based upon a composite figure of merit derived 
for each measurement system it is possible to relate the ex- 
pected accuracy of alignment to the expected accuracy of measure- 
ment. This includes the necessary conversion of linear measure- 
ments and magnitude measurements to their angular equivalents. 
The derivation of the composite figure of merit will be shown in 
what follows. 

2. 4 Error Analysis 

2.4.1 Angular Quantities 

In view of our basic definition of alignment, it is 
natural to speak in terms of error angles as an index of per- 
formance. All error sources contributing to the misalignment 
of a dependent reference frame may be described in terms of 
misalignments between appropriate pairs of the four basic coor- 
dinate frames of Section 2.2. Additional coordinate systems are 
used as analytical aids when necessary but are not central to the 
problem. 



The magnitude of any quantity properly described 
as a vector is invariant with respect to the coordinate transforma- 
tions necessary to characterize it in a particular reference coor- 
dinate system. The inputs to alignment measurement processes, 
as vectors, certainly have this attribute. If the errors resulting 
from a study of the measurement processes can be shown to have 
this attribute also, then to a very useful degree the alignment 
problem may be studied free of the geometry involved in specific 
hardware mechanizations. Because of the increased emphasis on 
inter -platform phenomena associated with the moving base align- 
ment problem, the choice of variables influences both the com- 
plexity involved in the analysis procedure and the insight afforded 
by the results. It is desirable that the variables chosen be appli- 
cable to all coordinate systems and thereby be independent of any 
particular mechanization. 

While angles cannot be described properly as vec- 
tor quantities, the simplest description of a misalignment angle 
nonetheless has invariant properties under coordinate transforma- 
tions. The simple description referred to reduces any defined 
difference in the angular orientation of two bodies to a single ro- 
tation about a specific axis. That this reduction may be achieved 

12 

is a consequence of Chasle's Theorem regarding the generalized 

1 5 

motion of a rigid body and has been discussed in detail in Ball's 
treatise on screw motion. 

In order to distinguish the single rotation, ® from 
its coordinatized representation or the series of distinct rotations 
employed to arrive at a final orientation, it will be called a "Whole 
Angle Rotation, " Although $ has the vector characteristics of 
magnitude and a directional axis, these properties are only ne- 
cessary to establish the whole angle rotation as a geometric entity 
in space whose representation is independent of particular coor- 
dinate systems. It will be helpful for the reader to visualize 
the alignment process as one of finding an axis about which one 
"turns" the dependent reference frame by a finite angle in order 
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NOTE: 

r * REFERENCE FRAME 
d * DEPENDENT FRAME 
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Figure 2-4 Pictorial Representation of the Whole Angle Rotation 

Required for Alignment 



to achieve angular coincidence with the reference coordinate set. 
This is schematically represented in Figure 2-4. The depen- 
dent reference frame, d, and the rectangular parallelopiped to 
which it is attached are rotated in space from their identical 
counterparts, the reference frame, r, and its associated 
parallelopiped. The alignment problem seeks only to achieve 
angular coincidence between the d and r frames. This may 
be done in several ways. 

In all cases of Euler type rotation, the following 
information must be recorded: 1) the body which is rotated; 2) 

the coordinate frame in which the rotations are specified; 3) the 
order of the three required rotations; and 4) the magnitude of 
each of the three rotations. On the other hand, if the whole an- 
gle representation of alignment is employed the following must 
be recorded: 1) the body which is rotated; 2) the coordinate frame 

in which the axis is specified; 3) the axis of rotation; and 4) the 
magnitude of the rotation, $. In terms of simple "bookkeeping", 
the Euler rotation requires keeping track of eight pieces of infor- 
mation; the whole angle rotation requires only five. 

The angle £>, in addition to being an easily visu- 
alized quantity, describes the magnitude of the misalignment by 
a single number. Until one requires a coordinatized representa- 
tion of this angld, the location of the rotation axis with respect 
to a particular reference coordinate set is unimportant. The mag- 
nitude <f> is, in fact, an upper bound on the magnitude of the angle 
<j)^ between any two corresponding misaligned axes. An original 
proof of this bound follows. The geometric definition of the </> . is 
depicted in Figure 2-5. <j> and are there seen to be the 

arguments of the principal direction cosines. 

That $ is an upper bound on the <j>. may be proved 
with the help of Figure 2-6. There the d and r coordinate systems 
are represented by sets of unit vectors, r. and d^. Beginning with 
the d and r frames coincident, the d frame has been rotated 
through the angle $ about the axis indicated to achieve the 


24 




orientation shown. Following, for example, the d^ and r^ axes, 
they are separated by the angle 0 ^ as a result of rotating the en- 
tire d-coordinate frame by $ relative to the r-frame. In general, 

= cos (r.. ’ d..) i = 1, 2, 3 2.4-1 

The rotation of d^ is seen most easily by decomposing r^ (co- 
inciding with the non - rotated d^) into a component perpendicu- 
lar to the axis of rotation, (r^)^, and a component parallel to 
the axis; (r^) (( . When the d-coordinates undergo the rotation 
$. (rj)^ rotates through the angle $ to the position indicated by 
(d^^while (r^)|| remains fixed. The displaced d^ axis is then 
constructed by the vector sum of (r^)j, anc * 

The geometry of Figure 2. 6 can be used directly 
to establish that bounding $ indeed sets an upper bound on 0^. The 
angles $and 0^ are each apex angles of two isosceles triangles 
which have the common base b. These triangles are drawn in 
the same plane in Figure 2. 7 for clarity. The triangles have 
sides [r^, H^]and L(r^)_, (d^jJ respectively. Since | rj = 
dj =1 and (r^)_J = I <1 we see immediately that 




2. 4-2 


where the equality holds when r^ is perpendicular to the axis of 
rotation. 

To obtain an explicit relationship between $ and 0, 
apply the law of sines to both triangles (Figure 2. 7) 


sin$ 


l<h>- 


sin a 


2. 4- 3a 


sin 0 


sin /3 


2. 4 -3b 
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Figure 2-7 Isosceles Triangles of Figure 2-6 Drawn in Same Plane 



Recognizing that 


sin 


sin /3 


1 *1 

sin *2 ( it - <^) = cos -g -1 2. 4-4b 


yields 


sin 




Lfjjj 

l f il 


sin 


$ 


2. 4-5 


If we define as the angle between the axis of the whole angle 
rotation and the coordinate axis of interest (in this case d^), then 


r l -L . 

. = sin 

' F 1 I 


2. 4-6 


Since the foregoing development holds independent of the coor- 
dinate considered. Equations 2.4-5 and 2.4-6 give the general 
result 




4> 


sm -g 1 = sin »|/ sin 


2. 4-7 


Figure 2. 8 shows <j> ^ as a function of $ for fixed values of ‘k. 
For very small <j > . and <£, the linear relationship 2.4-8 holds. 


<j)^ = (sin 40 $ 


2. 4-8 


$ as an exact function of the <j>^ is obtained by 
squaring Equation 2.4-7 and summing over i 
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L 

i= 1 
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Recognizing that ^ sin^ \\> ^ = 2 [derived from ^ cos^ vJj ^ 

i= 1 i= 1 


i =1 


Y" 2 2 

- 1 and ) (cos i|< . + sin ij; = 3] gives the desired relationship 

3 

2 $ i v . 2 h 

^ = Y 2 sin 2~ 


sin 


2.4-9 


i= 1 


When the small angle approximations are valid for $ and Equa- 
tion 2.4-9 or Equation 2.4-8 yields 





2. 4-10 


The relationship 2. 4-10 will be derived again further on in this 
section, without resort to the geometry of Figure 2. 6. From the 
result 2. 4-10 one may straightforwardly deduce Equation 2.4-2 
but not the explicit relation 2. 4-7. However, the forthcoming de- 
rivation demonstrates some important properties afforded by the 
use of the error angles </k. 

First the relationship of the magnitude of the whole 
angle rotation to several coordinatized representations commonly 
used in the analysis of inertial navigation systems will be computed. 
This is most easily accomplished by using matrix representations 
for the rotations. A dependent coordinate set, d, and a reference 
coordinate set, r, may be related by 



2. 4-11 


where R is a three -by-three orthogonal matrix whose elements 
are functions of the angular changes required to rotate the d- 
frame into coincidence with the r-frame. The reduction of R to 
a form which relates the r and d frames by an equivalent whole 




us- 


angle rotation proceeds by a similarity transformation 
ing a derived matrix S and its inverse S ^ 

R' = S R S" 1 

This equivalent rotation can be made of the form 

cos $ sin# 0 


12,16 
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-sin # 
0 


cos # 
0 


0 
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2. 4-12 


2. 4-13 


permitting evaluation of the magnitude of the whole angle rotation, 
#, from knowledge of the trace of R' 

t r R' = 1 + 2 cos # 2. 4-14 


Because the trace of a matrix is invariant under a similarity 
transformation, tr R = tr R' and therefore 

t r R = 1 + 2 cos $ 2. 4-15 


By 2. 4-15, the magnitude of the whole angle rotation,# , is evalu- 
ated from the trace of the matrix describing the desired rotation. 
This expression will be utilized frequently. 


Consider now the set of Euler angles corresponding 
to the rotations of a typical three- gimbal inertial system. A par- 
ticular set of rotation angles, 0^, and an order of rotation are 
defined in Figure 2. 9. Note that this set is the same as the tra- 
ditional yaw, pitch and roll of aircraft axes if dg, dg and d^ cor- 
respond to the aircraft's vertical, transverse and longitudinal 
axes respectively. The R matrix of Equation 2. 4-11 is, for this 
order of rotation, the matrix product 



where the abbreviations employed are 


sin 0. and C. = cos 0.. 
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SEQUENCE SHOWN IS 


Figure 2-9 An Euler Angle Transformation Corresponding 
to A Set of Three Gimbals 



Figure 2-10 Permutations of the Order of Rotations for 
Eulerian Transformation of Figure 2-9 



For the same values of the angles 9 ^, 6 ^ and 9 
there are five orders of rotation different from that described 
by Equation 2. 5-16 and Figure 2. 9. The cyclical permutations 
of Figure 2. 10 portray the entire selection of rotation sequences 
(Equation 2.4-16 is included in the cycle of Figure 2. 10 a). 
Although the matrix products LMN, MNL and NLM are unequal 
in general, the traces of these products are equal. The trace 
of the triple matrix product LMN may be written in terms of 
matrix elements as 


tr LMN ^ 

i j k 

that of MNL as 


t r MNL = 

■ m 

i j k 

and so forth, proving that 


III 
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1 . . m n. . 
i] jk ki 


m. . n., 1, . 
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2. 4-17 


1. . m n. . 
]k ki 


2.4-18 


tr LMN = tr MNL = tr NLM 2.4-19 
It may similarly be shown that 

tr LNM = tr NML = tr MLN 


■l 


i 3 k 


l.. n., m. , 
13 3 k ki 


2. 4-20 


and therefore 


tr LNM* tr LMN 


2. 4-21 


By Equation 2. 4-15, any matrices with equal traces 
describe equal magnitude rotations, We have thus shown that 
the magnitude of the whole angle rotation for any permutation 
within the same cycle of Figure 2. 10 is identically the same. 
Only the sense of the permutation affects the magnitude of <$, not 
the explicit order of rotation. This result is independent of the 
magnitudes of the 9 ^. 



The trace corresponding to the rotations of Figure 
2. 10a is (from Equation 2. 4-16) 

trR = trLMN = cos 0 ^ cos 9^ + cos cos ^3 

+ cos 0 g cos - sin 0^ sin 6 ^ sin 0^ 2. 4-22 

while that corresponding to the rotations of Figure 2. 10b is 
trR = trLNM = cos 0^^ cos 9^ + cos 0g cos 0g 

+ cos 0g cos 0 1 + sin 0^ sin 0 ^ sin 0g 2. 4-23 

Equations 2. 4-22 and 2. 4-23 differ only by the sign of the 
sin0^ sin 0g sin 0^ term. For the case of all 0. small enough 
to neglect terms of order three and higher, the traces become 
equal and reduce to 

tr LMNr tr LNM r 3 - 0j 2 - 0 £ - 0g 2 2. 4-24 


which with Equation 2. 4-15 gives 



+ 6 . 


+ 9, 


2.4-25 


The small angle result 2. 4-25 is completely inde- 
pendent of the order of the Euler rotations about the d^, dg and 
dg axes. This agrees with the well known result that the angles 

commute when dealing with "small" rotations. Markey and 
16 

Hovorka distinguish this situation by considering the "correc- 
tion angles" C ’ C , C in their work. From 2. 4-25 we imme- 
® x y z 

diately write 

<J>=/c 2 + C 2 +C 2 • 2. 4-26 

v x y z 


For errors specified by an azimuth error angle, 0 , 

a 


and a vertical error angle 0 , 
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2 0 2 + 0 2 
v a 


2. 4-27 
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While the trace of the rotation matrix has been 
shown to be relatively independent of the order chosen for a given 
set of Euler angles, it is quite dependent on the combination of 
Euler angles selected. Figure 2. 11 results from the following se- 
quence of rotations: Beginning with the r and d frames coincident, 

rotate the d-frame first through 0 3 about d^, then through 0 ^ about 
the displaced d^ axis, and finally through 9 ^ about the displaced d^ 
axis. The R matrix for this rotation has the form 

R = N (0 1 ) L (0 2 ) N (0g) 2.4-28 

where the L and N matrices are those of Equation 2. 4-16 with the 
appropriate substitution of trigonometric arguments indicated in 
Equation 2. 4-28. The trace of this rotation happens to be the 
exact trace for every permutation of rotational order. It is given 
by 


t r R = cos 0^ cos 0g + cos 0 2 + cos 0^ cos 0 2 cos 0g 
- sin 0 ^ sin 9^ cos 9^ - sin 0^ sin 0 3 2. 4-29 

which for small angles results in 

$ = yV + 0 2 2 + 0 3 2 '+ 20 1 0 3 2. 4-30 

Relations between $ and other Euler rotations can be evaluated 
similarly. 

The previous result 2. 4-10 can now be derived very 
easily. The arguments, of the principal direction cosines 
were defined by Equation 2. 4-1 as 

<f> i = cos" 1 (r. . d.) i = 1, 2, 3 2.4-1 

Because the dot products m . d^ are precisely the diagonal terms 
of the R matrix relating the d and r frames, the trace of the R 
matrix is written directly (equation 2. 4-31) 

t r R = cos </>^ + cos </> 2 + cos (j > 3 2. 4-31 
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The expression 2. 4-31 is exact and has been obtained without 
regard for the rotations or order of rotations involved. Any set 
of <j>^ describe a unique magnitude misalignment whose whole 
angle equivalent is -given by 

cos <£, + cos + cos <j > o -1 
cos $ = 2 ^ 2.4-32 

which reduces to Equation 2. 4-10 for the small angle case. 

+ *2 + *2 2.4-10 

In addition to the ease with which Equation 2. 4-31 
is written when dealing with the <f > this choice of error angles 
is such that the results 2. 4-32 and 2. 4-10 hold exactly even when 
the coordinate systems under consideration are not orthogonal. 

This result is useful in relating misalignments to measurement 
capabilities, particularly when a non- orthogonal reference frame 
is established with two of its axes along the directions of the two 
non-collinear vectors measured for alignment. Since measure- 
ment capabilities are determined as the ability to locate the direction 
of a vector, the individual misalignment evaluations remain un- 
changed as the angle between the two vectors changes. In addi- 
tion, if a .d-frame is chosen coincident with instrument axes, 
the form of the error analysis is relatively insensitive to input 
axis misalignment when carried out in terms of the $ 

To prove the independence of $from orthogonality 
of coordinate frames, consider the rotation of a non- orthogonal 
frame (n) to an orientation (n 1 ) such that 

sjc 

n' = R n 2.4-33 

where the asterisk distinguishes R as a non- orthogonal matrix. 

This can be reduced to the equivalent rotation of an orthogonal 
coordinate set (o) defined by 

o = A n 

o' = An' 2. 4-34 


37 



The rotation between triads is 


o' = (A R* A -1 ) o 2. 4-35 

(A is guaranteed non- singular provided no two axes of n are 
collinear. ) The similarity transformation and trace arguments 
of Equations 2. 4-12 through 2. 4-15, when applied to Equation 
2. 4-35 yield 

t r S(AR' A *) S 1 = t r AR* A 1 = t r R' = 1+2 cos $ 

2.4-36 

The trace of R is most easily expressed by choosing the angles 
between axes, </k, as the independent variables of R . Then 

3 

# ST 

tr R = ) cos 2.4-37 

1 

where the definition of Equation 2. 4-1 must be modified to read 

4>. = cos" 1 [m . EJ] 2.4-38 

The distinction is simply that the <£.. must be the angles between 
corresponding axes of two identical coordinate systems. If the 
d and r coordinate sets are defined alike, then Equation 2. 4-37 
is precisely Equation 2. 4-31 and Equation 2. 4-38 is the same as 
Equation 2.4-1. 

To quantify this distinction, consider two coordinate sets 
(n) and (n^) which are both non -orthogonal as well as non- identical. 
They are related by 

n = B 2. 4-39 

which establishes a coincident reference orientation of the n- 
frame and an n- equivalent n^ -frame. Keeping the n-frame fixed, 
rotate the n^ -frame to a new orientation n^' 

n^' = R' n 2.4-40 
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where R is now the rotation matrix between non- orthogonal, non- 
identical frames. Using 2. 4-39 

n 1 ' = R* B n x 2.4-41 


so that 


1 + 2 cos * = t r (R* B) 2. 4-42 

The trace of R B is a simple sum of the principal direction co- 
sines only if B is diagonal, corresponding to the n and n^ frame 
identical. 


The 0. will henceforth denote angles 7 between 
axes of identical (but not necessarily orthogonal) coordinate frames 
to permit evaluation of 4 from 


COS 0. + COS <j) 0 + cos - 1 
cos $= 2 — — 2.4-32 


or 



as applicable. 


2. 4-10 


The foregoing analysis of angular quantities 
which conveniently describe coordinate frame transformations con- 
cludes the first portion of the error analysis for this unified 
theory of alignment. As depicted in Figure 2-1, every align- 
ment procedure involves the transfer of a reference frame's orien- 
tation to a dependent reference frame. The quantity of interest 
for an error analysis is the misalignment between coordinate 
frames resulting from this transfer. Because this alignment is 
mathematically indistinguishable from a general rotation, the whole 
angle description has been developed in detail as the simplest 
description for misalignments. The whole angle description is 
especially convenient for dealing with the series of transfers in- 
voled in certain alignment techniques. 
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2. 4. 2 System Measurement Capability 

The actual transfer of coordinate frame orientation 
is accomplished by means of measurement processes. Errors 
made in these measurement processes are the fundamental 
source of the misalignments discussed in the previous section. 

In this section we establish the relationship between measure- 
ment errors and the resulting alignment errors. 

Vectors are the quantities measured in order to 
transfer alignment. As stated in Section 2. 3. 1, the primary 
object of measurements made for alignment purposes is to deter- 
mine the directions of these vector inputs. Although two non- 
collinear vectors are required to establish alignment unambigu- 
ously (Section 2. 3. 2), the measurement processes associated 
with these vectors are sufficiently independent of one another 
that the measurement of each vector may be considered separ- 
ately. As a consequence of these facts, the basic measurement 
problem may be stated as a simple question: "With what accur- 

acy can the measurements made with sensors integral to a par- 
ticular inertial system indicate the direction of a vector?" 

The phrasing of this question is such that the answer 
must be given in a particularly convenient form, namely as a 
single angular magnitude. This choice of an angular quantity as 
the interface between the measurement portion of the problem 
and the resultant alignment capability is a purposeful choice. 
While an angle is not necessarily the most natural means for spe- 
cifying the measurement capability of all types of sensors, it 
is indeed the most mathematically convenient for alignment dis- 
cussions. Additionally, it serves to highlight, a logical cause 
and effect pattern. The "effect" is always angular (misalignment) 
in nature while the "causal" factor may well be a linear measure- 
ment error. This latter situation is to be considered in detail 
shortly. Finally, the angular tolerance within which an instru- 
ment or set of instruments can indicate the direction of an input 



vector serves as a very natural "figure of merit" by which to 
compare various systems. The crux of any alignment ^technique 
is how well the direction of each of the required vectors can be 
indicated. 

The problem of vector direction indication has an 
exceedingly simple, yet important, pictorial representation. 

This is given by Figure 2-12. There, ly is the input vector nor- 
malized to unit magnitude, a is the angular uncertainty with 
which the measurement process indicates the direction of the 
unit vector, ly. In general, the plane defined by ly and a unit 
vector along the direction of the misindication will have no pre- 
ferred orientation about ly. This is an extension of the fact 
that even with perfect measurements, angular rotations about 
a vector cannot be detected by measurements of that vector alone. 
For the case of a non- zero measurement error, a is thus geome- 
trically interpreted as the half- cone angle of the cone whose gene- 
ratrix sweeps out the locus of possible incorrect indicated locations 
of the measured vector. The axis of the cone is along the true vec- 
tor direction. The apex of the cone is at the center of measurement. 
The location of a particular incorrect indication is singled out by 
the aximuth angle, rj. r\ is defined from an arbitrary (but fixed) 
reference and has any value on the interval 0 - 2n, all of which 
are equally likely. 

a is properly a random variable, not a determinis- 
tic parameter. This means that a actually describes a statisti- 
cally distributed family of coaxial cones, rather than a single 
cone. In most cases, the distribution of a is functionally related 
to the distribution of several other random variables. However, 
it is a single parameter which has a significant geometric meaning 
and which may be considered as a fixed quantity for purely geo- 
metric considerations. 

The uniform distribution of 17 is independent of the 
distribution of a . Yet it is recognized that the statistics of a may 


vary according to the relative vector-sensor orientation. Sta- 
tistical variations in both the magnitude and direction of a are 
logical. The magnitude variation of a is accounted for in this 
work. However, the fact that the misindication may be developed 
in a preferred direction relative to the sensor orientation is 
ignored here. Remember that our viewpoint is vector oriented, 
not system axis oriented. The preferred direction problem 
may be viewed as a second order refinement to be considered in 
the context of a detailed system analysis. 

Before proceeding to relate a to specific types of 
measurement errors, its role as a figure of merit requires 
clarification. Because we consider only the errors introduced 
by uncertainties in the measurement process, a actually bounds 
the ultimate accuracy of alignment which may be achieved with 
a particular set of sensors. It represents the best one might 
hope to achieve given a basic uncertainty in knowledge of the 
sensors. Unaccounted for is any slight degradation in accuracy 
introduced by processing of the basic measurement within the 
system. However, this bound emphasizes the basic limitations 
on alignment accuracy implied by varying measurement capabili- 
ties. 

2. 4. 2. 1 Basic Sensor Er rors and Their Equivalent 

Misalignments 

Our discussion of measurements concerns 
only "black-box" single-degree-of-freedom sensors. This re- 
striction is merely a conceptual aid which does not preclude the 
use of multi-degree-of-freedom sensors in an actual system. 

Two- or three- degree-of-freedom sensors are mathematically 
equivalent to two or three single-degree-of-freedom sensors 
appropriately mounted on a rigid structure. 

There will follow shortly an examination 
of the ability of various groupings of these single-degree-of- 
freedom sensors to indicate the direction of a vector input. 


First, however, we undertake a more abstract (but instructive) 
examination of the measurement problem by briefly considering 
a hypothetical instrument which aligns itself along a vector in- 
put (a "vector seeker", if you will). Every realistic sensor is 
characterized by two distinct quantities: 1) a sensitive axis 
(fixed relative to the sensor's case); and 2) a certain transducer 
behavior with respect to this axis. We will so characterize the 
"vector seeker". 


Due to errors in the construction and 
mounting of an instrument, its sensitive axis is typically not 
located with perfect precision. We call this imprecision, "sen- 
sitive axis uncertainty" and denote it by the angle a g^; so too 
for the vector seeker. Imperfect transducer behavior of an in- 
strument also results in an angular "indication error, a j,- This 
we denote transducer uncertainty. Again, this nomenclature is 
carried over for the vector seeker. For clarification of the mea- 
ning of the basic instrument uncertainties, the reader is referred 
to the more classical example of Appendix A. 


The virtue of developing the vector seeker 
concept is that it reduces various mechanizations of the same 
basic problem to a common denominator. Our ability to evaluate 
and compare different sensor systems is thereby greatly enhanced. 
The development of the total indication uncertainty, a, of Figure 
2-12 from ct and a now proceeds for the vector seeker. 

The "vector seeker" uncertainties equivalent to appropriate ar- 
rangements of single-degree-of-freedom sensors are calculated 
later in this section. 


A vector seeder which is a perfect trans- 
ducer indicates a signal null when its sensitive axis is collinear 
with a vector input. The geometry of such a situation is depicted 
in Figure 2- 13a. As long as an observer knows the precise lo- 
cation of the sensitive axis, the vector's direction will be found 
without error. On the other hand, in spite of the fact that the 
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"VECTOR 

SENSITIVE 

SEEKER" 



VECTOR INPUT; 
INDICATED INPUT 


(a) PERFECT SENSOR 



(b) SENSITIVE AXIS UNCERTAINTY 
(LOCUS OF CONSTANT PROBABILITY DENSITY) 



(c) TRANSDUCER UNCERTAINTY 
(LOCUS OF CONSTANT PROBABILITY DENSITY) 


Figure 2-13 The Geometry of Vector Direction Indication 





(d) TRANSDUCER AND SENSITIVE AXIS UNCERTAINTY COMBINED 


A 



(e) SPHERICAL TRIANGLE FOR RELATING a, a T/ & a $A 


Figure 2-13, continued. 
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true instrument sensitive axis lies along the vector, any sensitive 
axis uncertainty (u g^ of Figure 2- 13b) is directly reflected in 
an equal magnitude direction indication uncertainty. The con- 
stant probability density locus for misindications due to uncertainty 
in the location of a sensor's sensitive axis is a right circular 
cone whose axis is the true sensitive axis and whose half-cone 
angle is a g^. Because the magnitude of a g^ is dependent 
merely on a defined quantity, it has no influence on the physi cs 
of the measurement itself. It is a random bias angle which can- 
not be calibrated out. 

Figure 2-1 3c illustrates the situation where 
a g^ = 0 but the transducer behavior of the "vector seeker" is 
imperfect. These errors affect the dynamics of the measurement 
process, resulting in the "seeker's" true sensitive axis coming 
to rest a ,p from the vector's direction when a signal null is 
reached. The constant probability density locus for misindications 
due to imperfect transducer behavior is a cone whose axis is along 
the vector input, whose half- cone angle is a -p, and whose gene- 
ratrix is the true sensitive axis of the sensor. 


Combining the effects when both a g^ and 
a r-p are non-zero results in Figure 2-1 3d. The total indication 
error, a , is the a of Figure 2-12. It is measured between the 
generatrix of the cone of half-angle Ug^ and the direction of the 
vector input. The axis of the cone of half-angle is the gene- 
ratrix of the cone of half-angle u-p. Note, therefore, that the 
surface of the cone of half-angle Ug^ is the locus of misindicated 
directions only for the illustrated location of the true sensitive 
axis. The complete locus consists of the surfaces of all cones of 
half-angle Ug^ drawn about the infinity of allowable positions of 
the generatrix of the cone of half -angle a ,-p. 

In spite of the seeming complexity of this 
locus, the evaluation of a is straightforward. An exact relation- 
ship between the contribution of transducer uncertainty, a^,, and 


46 



the contribution of sensitive axis uncertainty, to the total 

indication error, a, is obtained by analyzing the geometry of 
Figure 2-1 3d on a unit sphere. The center of the sphere is 
coincident with the center of measurement.. The vertices of the 
spherical triangle (Figure 2-1 3e> of sides a^, and a are found 

by the intersection of the following three radial directions with 
the surface of the sphere: the vector input direction (vertex A); 
the direction of the true sensitive axis (vertex B); and the di- 
rection of the assumed sensitive axis (vertex C). Defining the 
angle p between or^, and allows writing the law of cosines for 
spherical triangles 


cos a = cos cos <*g^ + s i n sin ^SA cos ^ 


2. 4-43 

For small values of a ^ and a g^. Equation 2. 4-43 becomes 

2 


2 _ 2 
a a T + a SA 


- 2 o' T cos p 


2. 4-44 


17 may be written in terms of the azimuth angles 17 ^ and p 2 the 
cones of half angle cv,p and <*g^, respectively 


f] = * - (r ) 1 - p 2 ) 


2. 4-45 


Substituting Equation 2. 4-45 into Equation 2. 4-44 yields 
2 2 2 

a = ar T + tfg^ + 2tt q' a g^ cos 2 * 4 "46 

p^ and p 2 are statistically independent. Each is uniformly dis- 
tributed on the interval 0 - 2tt. 


Equation 2. 4-46 is useful as written only 
when exact values of p^ and p 2 are known. This corresponds to 
knowing the relative orientation of the true and assumed sensitive 
axes and is in violation of our basic error model. Meaningful 
quantities which may be readily ev alua ted from Equation 2. 4-46 
are the mean (a) and mean square ( a ) values of a . 

The mathematics for this will be found in 
Appendix B, with the results 
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2 . 4-47 


a ~ \ (« T + a g^) £ (k 2 ) 

2 2 , 2 „ , , 0 
a = + a gA 2.4-48 

1 8 

where £ is a complete elliptic integral of argument 


2 

(ar T + “SA> 2 


Note that the results of 2. 4-47 and 2. 4-48 are purely geometric. 
They were obtained without knowledge of the statistical distribu- 
tions of and Ug^. The mean and mean square values of a as 
functions of the statistics of and Ug^ are obtainable from Equa- 
tions 2. 4-47 and 2. 4-48 by further evaluation of the expectations 

indicated there. While a cannot be evaluated without knowledge 

2 

of the distribution of and a g^, a is immediately given by 


, a 1 - + <*g A 2 2.4-49 

As stated at the outset of this section, the 
single-axis "vector-seeker" merely represents a system of from 
one to three single-degree-of-freedom sensors. In order to utilize 
this concept, we must derive the Q'g^ an< ^ related to mea- 

ningful groups of sensors. 

We first set u . to zero in order to inves- 

O-A. 

tigate the relationship between and various types of transducer 
inaccuracy. In certain cases, the basic output of the sensor is 
proportional to an angular quantity, one of whose rays is the 
measured vector. Optical instruments, such as star trackers and 
auto- collimators, fall into this category. In other cases, the sen- 
sor's output signal is proportional to the magnitude of the vector 
and this signal must be further processed to establish the vector's 
direction. Inertial sensors, such as gyroscopes and accelero- 
meters, are in this category. In the former situation, error 
quantities are normally specified as angular deviations from a 
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nominal input direction. These specifications would determine 
ctrj, directly. In the latter situation, however, error quantities 
are related to magnitude uncertainties and these, in turn, must 
be converted to equivalent angular uncertainties. After conver- 
sion, the magnitude measurement errors determine a ^ 

There are two distinct procedures for 
using magnitude sensitive instruments to determine a vector's 
direction. The first is to determine the direction directly 
from magnitude information. The second is to determine the 
null plane of the vector and infer the vector's direction from 
this. In the following, consider at first the measurements made 
by a single instrument and assume a priori knowledge of the mag- 
nitude of the vector or the significance of its null plane. 

The cases of estimating the magnitude 
of a vector and estimating the null plane of the vector must be 
considered separately. Associated with these measurements 
are the magnitude estimation error quantities e m and e n respec- 
tively. e is the error in the indicated measurement when the 
instrument's input axis is along the input vector and the input 

magnitude has been normalized to unity, e is the indicated out- 

■ n 

put when the input axis is in the null plane (the ideal output is zero) 
normalized against the magnitude of the vector quantity whose 
null is sought. 

Figure 2-14 illustrates the case of attempting 
to measure the orientation of a normalized vector, ly, by mea- 
suring its magnitude. If the instrument is high reading (indicated 
input exceeds actual input by e ), it will indicate a unit magni- 
tude when, in fact, there is only a projection of (1 - e ) along its 
input axis (Figure 2-1 3a). This corresponds to the angular error 
Oj. indicated. If, on the other hand, the instrument is low reading, 
the relationship of Figure 2-1 3b holds. When the instrument is 
sensing an actual input magnitude (1+ e ) it indicates that, it is 
along the unit magnitude vector Ty and the desired direction is 
removed from the indicated direction by 



PROJECTION OF 
K\ LENGTH |1-« M ) 


DESIRED 

DIRECTION 


INDICATED 

DIRECTION 



DIRECTION DIRECTION 


(o) 


lb) 


HIGH READING INSTRUMENT 


LOW READING INSTRUMENT 


Figure 2-14 Vector Magnitude Measurement Errors 



Figure 2-15 Vector Null Plane Measurement Errors 
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For the high reading case 


cos a = 1 - e 2. 4-50a 

a m 

while for the low reading case 

sec a, ; 1+e 2. 4- 50b 

b m 

Using the series expansion for the cosine and ignoring terms of 
order greater than two in a yields, respectively 

2. 4-51a 

2. 4-51b 
direction of ly 


“a = / 


2 e 


m 


and 



m 


so that when e 1 the error in indicating the 

m & 

by measuring its magnitude becomes 


a 

m 


a ~ a 


a = / 2 e 2. 4-52 

m J m 

The measurement of a vector's null plane 
is shown in Figure 2-15. The desired normal direction is ly, 
but due to the existence of e n normal to the unit vector ly, the 
indicated normal direction is rotated from the desired normal 
direction by an angle o^. The null planes are therefore also <* n 
apart. Geometrically, tan = e n , which for small e n reduces 
to 

a = e 2.4-53 

n n 

The angles and t* n are the a-^ for each case. They represent 
the necessary conversion of magnitude errors to angular equiva- 
lents . 



Equations 2. 4-52 and 2. 4-53, although 
derived for a restricted situation, contain an interesting funda- 
mental result. For a fixed measurement capability (e = e 
but e <2) greater alignment accuracies can be achieved by me- 
chanizations operating on a nulling principle (Equation 2. 4-53) 
than by those operating on a magnitude measurement principle 

(Equation 2. 4-52). For typical values of e in the range of 10 ^ 

— 6 

to 10 , the alignment capabilities vary by one-half an order of 

magnitude to over three orders of magnitude. In many cases, 
e n < e , widening the gap even further. 

Unfortunately, the results 2. 4-52 and 2. 4-53 
present a picture which is oversimplified on two counts. First 
of all, a single instrument will rarely suffice for either magnitude 
measurements or null-seeking measurements. It is more com- 
mon in inertial systems to provide sufficient mutually orthogonal 
sensors to span the defined measurement space. Thus for null- 
seeking measurements, a minimum of two instruments is required 
to define the null plane. Except for special situations, magnitude 
measurements require three sensors. We must, therefore, de- 
termine the resultant direction indicating capability of appropriate 
groups of sensors rather than just that of single instruments. This 
will be calculated shortly. 

The second oversimplification is in our 
error model for the inertial sensors. The quantities and 
imply that an instrument is defined by a fixed percentage er- 
ror. This is not so for precision inertial sensors. Rather, the 
indicated output (V ^) of a device is more commonly expressed 
as a weighted power series of the input (V) along its sensitive 
axis plus terms dependent on the sensor's environment (thermal, 
field, etc. ) and terms related to inputs normal to its sensitive 
axis (See Equation 2. 4-54). 
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v 

ind 


+ a, V + 


a 2 v2 + 


+ fen of environment + cross axis coupling 

2. 4-54 

There is a popular nomenclature associated with the coefficients 
a^ of Equation 2. 4-54 which depends upon the instrument under 
consideration. Hence a^ represents gyroscope drift equally as 
well as accelerometer bias (See Appendix A. ) These distinc- 
tions are relatively unimportant in this paper as numerical 
alignment results are expressed in terms of non-dimensional 
quantities. 

For convenience, all derivations will 
assume a normalized instrument whose ideal transfer from in- 
put to output is unity. Thus a perfect instrument is described 
by Equation 2. 4-54 with a^ = 1 and all other coefficients zero. 

In this work, imperfect sensors will be described by the trun- 
cated representation 2. 4-55 rather 2. 4-54. 


V 


ind 


+ a t V 


2. 4-55 


The term "bias " will be synonymous with a^, "sensitivity" with 
a^. The uncertainty, .. or error, in a^ is defined as Aa^ = (a - 1). 
Because the error coefficients a^ and Aa^ of an instrument are 
generally compensated to some degree in a system application, 
we take a^ and Aa^ to represent the uncompensated uncertainty 
existing in the sensor at the time of measurement. They are the 
residual error between the actual value of the coefficient and the 
best estimate of the coefficient available to the system at the time 
of alignment. 

2.4. 2.2 System Magnitu de Measurements 

Let us now reconsider the directional in- 
dication error associated with magnitude measurements of a 
vector V (IV | = V). Although we are most interested in 

measurements made with a triad of instruments, it is instructive 


53 



to rewrite the single instrument result (Equation 2. 4-52) in terms 
of the description 2. 4-55. When indicating the direction of a vec- 
tor with a single instrument, the indicated vector magnitude (V\ 
is constrained to equal the true vector magnitude (V). In order to 
satisfy this constraint a directional error of 

-J2 ( a o/y + 2.4-56 

must exist. This error is physically realizable only if 
( a o/v> + ^ a l closed loop operation against an a priori mag- 
nitude is impossible to achieve with a low reading instrument. 

This is simply because the output of a low reading instrument can 
never be as large as the known magnitude of the measured vector. 
Only an open loop search for a maximum indicated magnitude is 
physically realizable on an unrestricted basis (It happens also to 
be free of transducer error in the steady state). Regardless of 
the restricted applicability of the expression 2. 4-56, it demon- 
strates that the angular uncertainty resulting from locating a vec- 
tor with a single instrument is an explicit function of V, a^ and 
Aa^. Contrast this to the implicit form of Equation 2. 4-52. 

The explicit form shows that (and thus the directional indi- 
cating capability) cannot be divorced from the magnitude of the 
input to the sensor and cannot therefore be generally described 
as a constant percentage. In fixed base alignment techniques, 
the viewpoint of the system designer is to take the magnitude of 
the measured vector (e. g. , | g | or |u ^ | ) as a fixed quantity and 

adjust his alignment capability, primarily through the choice 

of instrument parameters. In this context, once his choice of 
instrument is made e may be properly thought of as a fixed 
percentage error. By contrast, the designer faced with a moving 
base alignment problem must often consider the magnitude, V, 
of the vector as another variable over which he exercises con- 
siderable choice. In fact, by specifying the carrying vehicle's 
motion during alignment he chooses the dimensions, direction 
and magnitude of V. His choice is bounded only by the opera- 
tional restrictions on movements of the carrying vehicle. The 



accuracy achieved at any given time depends on just how 
closely V approaches a desired nominal value. Because V may 
not be considered fixed in moving base alignment, e m may not 
properly be considered a fixed percentage. 

We proceed now with the more realistic 
problem of indicating the direction of a vector, V, with an ortho- 
gonal triad of instruments characterized by 

a^ = Bias uncertainty of i *'* 1 instrument. 

"t It 

a^ = Scale factor of 1 n instrument. 

Aa^ = Scale factor uncertainty of i th instrument. 

- "If 1 

V is described in instrument coordinates by 

V = V (cos i + cos i[> g j + cos ijj 3 k) 2. 4-57 

while the vector indicated by the instruments, is described 

in the same coordinates by 

hnd = Cta 01 + an Vcos + pT 
+ (a 02 + a 12 V ” oos + 2 > j 

+ (a Q3 + a 13 V cos i|j g ) k] 2.4-58 

(See Figure 2-16) Equation 2. 4-58 may also be written 

hnd = V + ? 2. 4-59 

where 

e = + Aa.^ V cos ip i 

+ (a 02 + Aa i2 VCOS4/ 2 ) T 

+ ^ a 03 f Aa 13 ^ COS ^ 3^ ^3 2. 4-60 
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Figure 2-16 Vector Direction Indication With A Triad of Sensors 





It is reasonable to expect that on the average any worthwhile 
system of sensors will indicate V perfectly. That is, the mathe- 
matical expectation of Vj nd is V 

E <V ind ) = V 2.4-61 

and consequently 

E ( 7 ) = 0 2.4-62 

Equation 2. 4-62 may be guaranteed independent of the magni- 
tude or direction of V provided 

E (a n .) = E(Aa..) = 0 2.4-63 

Ui li 

The techniques for establishing a zero mean-value for sensor 

coefficients (via instrument design and calibration) are suffi- 

19 

ciently well established that this paper takes the statement of 
Equation 2.4-63 as a basic assumption rather than as a require- 
ment. 

The reader is cautioned that Equation 2. 4-61 

must be interpreted literally. It is a vector equality stating that 

the vector V^ nd equals the vector V on the average. This does 

not imply that, for instance, the magnitude of Vj d equals the 

magnitude of V on the average. Nor does it imply that the average 

direction of V^ nd is that of V. In spite of our intuition that the average 

measured magnitude and direction are correct it merely states 

that the average total measurement of V is correct. This point is 

20 

often missed in optimization studies which concentrate on tra- 
jectory endpoint miss- distance minimization. The distinction 
will be quantified shortly. 

For alignment, of course, the direction 
indicating capability of the triad is of greatest interest. An 
error in directional indication is expressed as the angle be- 
tween V and V\ Because the triad is insensitive to rotations 
about V, we ignore the triad-referenced direction in which ar.^, 
is developed and concentrate on evaluating the magnitude of ovp. 



This is an approximation which merits closer attention only 
when the platform rotation about V is precisely constrained. 


from 


The magnitude of a j, may be calculated 


COS QVj, 



( 

! 1 + 2 


111 + 111 ) 
V 2 V 2 J 


1/2 


2. 4-64 

which may be expanded to 



2. 4-65 

Note that although |e | « |V| , terms of e ' e may not be ignored 
in comparison with terms of e ’ V because orthogonality is suffi- 
cient to cause e * V = 0. 

For three instruments of equal quality with 
symmetrically distributed coefficients (e. g. , normal, uniform) 
and assuming (1) that each instrument is independent of every 
other instrument and (2) that the bias and scale factor uncertainty 
of any one instrument are also independent, it is a straightforward 
procedure to compute the mean- square value of a as 



2 


E(af T ) - 2 


V 



I 

- i 


4 


cos 



2. 4-66 


where terms of fourth- and higher-order smallness have been 

dropped from the result 2. 4-66. a and cr A are the standard 

a 0 Aa l 

deviations of instrument bias and scale factor respectively. 


Equation 2. 4-66 provides a good deal of 
insight to the direction indication problem. First of all, it is 
possible to minimize the mean- square value of by adjusting 
the orientation of the triad with respect to V. Referring to 


Figure 2-17, the function 



4 

cos 




is zero when V is 


parallel (or anti-parallel) to any one of the instruments' sensi- 
tive axes. Not only does this physical orientation correspond to 

null seeking with two instruments but also the mean-square 
2 o’ 2 

Etct'rp'j = 2 ( a^/V) is identically that derived for the null seeking 
case (See Equation 2. 4-78). This is interesting because the 
result appears to be independent of whether or not the instrument 
along V is operative or inoperative. Why, in fact, does the fac- 
tor 2 appear in Equation 2. 4-66 when three instruments are 
active? In answering this, remember that we are only consi- 
dering direction indication. The probability density for bias 
is spherically symmetric (for equal instruments), but only 
components of bias in the plane perpendicular to V contribute 
to u,p. The very important contribution of spherical symmetry 
is that the tangential component of bias always appears due to 
two instruments while the radial component appears ,due to the third. 
The radial component affects only the magnitude of "V\ ^ (See 
Equation 2. 4-67). 


From these considerations we draw the 
important conclusion: On the basis of minimum mean-square 
error, it is optimum to indicate the direction of a vector by 
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3 

Figure 2-17 The Surface Described by a Radius of Magnitude (1 - V cos^ .) 

Lj i 

i=l 


seeking its null plane with two instruments. This result strongly 

17 

supports the practical reasons usually advanced for placing 
the accelerometer input axes of a terrestrial navigator in local 
geographic coordinates. 

2 

The maximum value of E (a ) occurs when 
_ T 

V is symmetrically located with respect to the sensors, i. e. , in 

the center of the triad. It is relatively easy to show the impor- 
tance of orientation on the error As an example, evaluating 

Equation 2. 4-66 with constant a , cr A , and V demonstrates that 

a 0 a l 

merely by adjusting the location of the vector relative to the triad 

2 

there is at least a 10% reduction of E (a,j, ) available provided 




Aa, 


For 


4 s 


the reduction 


available by choice of orientation is less than 10 percent. In 
situations where the bias uncertainty predominates, therefore, 
the orientation of V with respect to a triad of sensors becomes 
less significant. 

We turn for a moment to examine the accuracy 
with which the preceding triad measures the magnitude of V. The 
mean- square magnitude measurement error may be calculated as 

3 


E C(l V. A - Ivl ) 2 ] 


ind 


+ V' 


! 2 V 

CT Aa 2j 

1 i=l 


cos 


+ 1 


2. 4-67 


As with Equation 2. 4-66, Equation 2. 4-67 ignores terms of 
fourth- and higher-order smallness. From Figure 2-18 it is 
evident that the minimum value of Equation 2. 4-67 occurs when 
V is in the center of the triad and the maximum occurs when V 
is along any one of the instrument axes. This is clearly opposite 
to the results for direction indication. See Table 2-1, The 
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applicable statement here is: On the basis of minimum mean- 
square error, it is optimum to measure the magnitude of a 
vector with a triad of instruments by placing the vector at the 
center of the triad. 


Vector Orientation 

Magnitude Error 

Direction Error 

Along any instrument 
axis 

Maximum 

Minimum 

Center of Triad 

Minimum 

Maximum 

1 


Table 2-1 

Demonstrating the Difference Between 
Magnitude and Direction Measurement 
Errors with a Triad of Instruments 

For any situation simultaneously requiring accurate measure- 
ments of the direction and magnitude of a vector, the choice of 
vector orientation with respect to the triad is obviously a com- 
promise. For example, the thrust vector is an attractive vector 
for alignment purposes but accurate knowledge of magnitude of 
this vector is simultaneously important for guidance. 

An interesting constraint exists between 
mean- square magnitude and direction errors. It derives from 
the fact that the mean- square value of the error vector e is in- 
dependent of the vector's orientation. See Equation 2. 4-68. 

E (e 2 ) = 3 a 2 + V 2 cr 2 * 2.4-68 

a 0 ^ a l 

Equations 2. 4-66 through 2. 4-68 may be combined to show 
E (? 2 ) = V 2 E (a T 2 ) + E [ ( |v. nd i - IV| ) 2 ] 

2. 4-69 
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For small a , M 2 - (Va ) 2 + (IV jnc| l - IV I ) 2 

e = total indication error vector 
a = _di rection indication error 
IV. | - | V f = magnitude indication error 


Figure 2-19 Relationship Between Magnitude and Direction 
Measurement Errors for a Given Triad 
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Figure 2-19 presents a graphical interpretation of this equation. 
It is accurate to the extent that the chord length subtended by 
the arc Vor numerically equals Vh,p. The available tradeoff 
between magnitude and direction errors is thus seen to be sim- 
ply described (Equation 2. 4-69) as a linear constraint. In a 
situation such as this, the discretion of the designer in choosing 
an appropriate compromise is extremely important. 

It is also interesting here to compare the 
quality of magnitude estimation attainable with a triad of instru- 
ments and that attainable with a single instrument whose input 
axis is perfectly aligned with the vector. On the basis of mean- 
square errors, this ratio is written as 


Magnitude Estimation Error of Triad 
Magnitude Estimation Error of Single Instrument 


1 + r 


cos 


*, • 


i=l 


1 + r‘ 


2. 4-70 


where 


CT 

a 0 

The upper and lower bounds of 2. 4-70 are plotted in Figure 2-20 
as a function of r. With any one instrument of the triad along 
the vector, the ratio is unity independent of the value of r. In 
other words, the triad is then no better than a single instrument. 
When the vector is in the center of the triad, however, there is 
an improvement with the triad of up to a factor of three (for 
predominant scale factor error). In between these extremes 
of no improvement and a factor of three improvement in magni- 
tude estimation, the parametric variation is dependent upon both 
the value of r and the relative location of the vector with respect 
to the triad. As with the direction indication problem, the 
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RATIO OF MEAN SQUARES 



MAX. ERROR; ANY INSTRUMENT 
ALONG VECTOR 


VALUES OF RATIO LIE IN 
SHADED REGION 


MIN. ERROR ; VECTOR IN 
CENTER OF TRIAD 


Figure 2-20 Comparing Magnitude Estimation Error of a Triad to 

That of a Single Instrument 



SINGLE AXIS SENSOR 


Figure 2-21 Relative to Vector Direction Indication by Means of 
Two Magnitude Sensitive Instruments 
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orientation of the triad is least significant when bias uncertainty 
predominates. 


One additional situation is relevant to 
direction indication via magnitude measurements. This entails 
indicating the direction of a vector by use of two orthogonal in- 
struments and is germane to azimuth alignment by a vector 
matching procedure (see Chapter III). Using the geometric 
definition of Figure 2-21 and assumptions similar to those em- 
ployed in the derivation of Equation 2. 4-66, we obtain the result 


a 0\ 2 


4 4 

1 - (cos vJj + sin 4 * ) 


This is rigorously correct for situations where V is out of the 
plane of the instruments provided V represents the magnitude 
projection onto the instrument plane. The error 2.4-71 is mini- 
mum when either instrument is along the vector fib = 0°, 90°, 
180°, 270°). In these orientations the scale factor uncertainty 


makes no contribution to the mean- square direction indication 
error. At = 45°, 135°, .... the maximum error of 
[<%/V > 2 + o 1 12 CTAa ?H occurs. There is at least a 10% reduc- 
tion of E (urp^ ) available by orientation provided (°aQ/y) — 


2. 12 <7 


Aa l‘ 


Results have been derived in this section 
concerning the errors resulting from indicating the direction of 
a vector with various groups of magnitude sensitive instruments. 
The mean- square results are combined with those of the next 
section concerning null-plane seeking with the same instruments 
and are presented in Table 2-3. 


2. 4. 2. 3 System Null-Seeking Measurements 

We turn now to the problem of direction 
indication by null- seeking measurements. The explicit object 
of null- seeking mechanizations is to determine the null-plane of 
a vector. Implicit in this determination is the knowledge that the 
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nonnal to the null-plane is parallel (or anti-parallel) to the di- 
rection of the vector. If it is necessary to determine the signed 
direction of the vector, the sign ambiguity inherent in null- 
seeking must be resolved by other measurements. 

The principle of null-seeking measure- 
ments is to control the orientation of a sensor's axis in order to 
achieve an indicated output of zero. For a perfect sensor, the 
null indication is achieved when the sensitive axis is normal to 
the vector input. This is independent of the vector's magnitude. 

A sensor with zero bias error but with a sensitivity error will 
likewise have zero output when its sensitive axis is in the null 
plane. Directional indication errors can arise only in the presence 
of an instrument bias. In order to obtain a net sensor output of 
zero when bias is present, the sensitive axis must be tipped out 
of the null plane an amount sufficient to cancel the bias term 
with a small projection of V along the sensitive axis. From Equa- 
tion 2. 4-55 the constraint is 

a^ + V sin = 0 2.4-72 

where is the angle between the sensor's sensitive axis and the 
normal projection of the sensitive axis on the null plane (i. e. , 
the complement of the angle from the sensitive axis to V). for 
small 



By comparison with 2. 4-53, the magnitude of the parameter e n is 
seen to be |e I = <an/v> . The fixed percentage error description 
n a l 

of inertial instruments is thus invalid for nulling mechanizations 
as well as magnitude measurements. 
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The complete indication of a vector's null 
plane requires two instruments, each of whose capability is de- 
scribed by Equation 2. 4-73. We calculate now the net transducer 
inaccuracy, < 2 ^, of Equations 2. 4-48 and 2. 4-49 for the complete 
null-seeking mechanization. In Figure 2-22 the instruments are 
distinguished by the subscripts a or b. 1 is a unit vector along 
the input vector direction. 1 & and 1^ are unit vectors along the 
sensitive axes, separated by the angle 0(0 = 90° nominally) in 
the plane of the sensitive axes. The indicated direction of ly 
is the direction of the cross product 1 x 1 K , shown by the unit 
vector a a anc * « b represent the a ^ of Equation 2. 4-73 for 

each of the instruments. 

is derived by writing 1 as components 

along T a , T b andT ab 



(1 V - 1 a ) + X b (1 V * 1 b ) + J ab (1 V ' 1 ab ) 


= 1 sin a + 1, sin u, + 1 cos 2. 4-74 

a a b b ab T 

Taking the inner product of ly with itself and solving for cos 
gives the result 


COS QTrp 



sin a, - 2 sin a sin a, cos 0 
b a b 

2. 4-76 


For small error angles this reduces to 

Qfrp + ar b ^ + 2 a & o? b cos 0 2. 4-77 

Both 2. 4-76 and 2. 4-77 emphasize the sign ambiguity inherent in 
null- seeking mechanizations. In fact, the a t of Equation 2. 4-73 
does not have to be carried as a signed quantity. Only its mag- 
nitude is important. 
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By 2. 4-77 a null-seeking mechanization 
with two orthogonal instruments of equal quality is found to have 
a resultant mean-square transducer error of 




2.4-78 


This expression recognizes that the error contribution of scale 

factor uncertainty in nulling mechanizations is of secondary im- 

2 

portance compared to the bias contribution. c*,p represents one 
of the statistical parameters of a "vector seeker" operating on 
the nulling principle. There is no spatial dependence of a ^ as 
with the triad of magnitude sensitive instruments, but this does 
independently check the results obtained for a triad (Equation 
2.4-66) when any instrument is along the vector. Since 2. 4-78 
was calculated on the basis of two instruments and 2. 4-66 on the 
basis of three instruments, it is indeed interesting that the re- 
sults are identical. 


It is instructive to investigate the relative 
importance of matching instrument quality versus maintaining pre- 
cise orthogonality. Fqr orthogonal but unequal instruments. 
Equation 2. 4-77 may be written as 


a T 


+ 1 a. 


2. 4-79 


a a % 

where n is chosen from either n = or n = — in order to 

O'. O' 

b a 

constrain n — 1. For equal but non -orthogonal instruments 


a T = J 2 + 2 cos 6 = (2 cos g ) 2. 4-80 

The points where instrument mismatch and non- orthogonality 

make equal contributions to the misalignment are obtained by 

equating 2. 4-79 and 2. 4-80. The percentage mismatch between 

instruments is given by 100 (n-1) while the percentage orthogo- 

Q 

nality error is given by 100 (1 - The curve of equal 


71 



% ERROR IN ORTHOGONALITY 
(1%= 0.9*) 


Figure 2-23 Curve Relating Orthogonality and 
Instrument Mismatch Errors whichCause Equal 
Misalignments 





misalignment contributions is plotted in Figure 2-23. This shows 
that on a percentage basis the total null- seeking capability is 
more dependent on maintaining orthogonality than on matching 
instruments. For example* an error caused by a 3% instrument 
mismatch is the same as that caused by a 2% orthogonality error. 

The results of this section concerning the 
errors resulting from seeking the null-plane of a vector with 
magnitude -sensitive instruments are tabulated in Section 2. 4. 2. 5* 
Table 2-4. 


2. 4. 2. 4 System Sensitive Axis Uncertainty 

In this section the equivalent sensitive axis 
uncertainty, Og^, for systems of sensors is calculated. The 
measurement principles considered correspond to those discussed 
in the preceding two sections. In keeping with the earlier assump- 
tion of equal quality sensors for each system, the instruments 
are all assumed to have the same sensitive axis uncertainty, a 

SSL 

Additionally, the sensitive axes of multiple sensor systems are 
assumed orthogonal. <*,j, is temporarily set to zero-. 

Both magnitude and nullipg measurements 
with single instruments are trivial cases for which the equivalent, 
uncertainty and instrument uncertainty are identical (cv OA = a ). 

uA SSL 

Two sensors employed in a nulling me- 
chanization, however, present a much more difficult situation. 

With multiple sensor operation, we revert once again to the 
"vector seeker" concept in order to define the equivalent < 2 g^ 
for the system. The angle between the assumed and actual sen- 
sitive axes of a system's equivalent "vector seeker" is the desired 
angle c*g^. For the nulling system, therefore, is the angle 

between the normal to the plane containing the actual sensitive 
axes and the normal to the plane defined by the assumed sensitive 
axes. Figure 2-24 a shows this situation. The unit vectors i 
and j lie along the instruments' actual sensitive axes while the 
unit vectors 1 & and 1^ lie along the assumed sensitive axes 
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(coinciding with the "scribed line"). The axes are related by 
the constraints 


1 . i = cos a 
sl sa. 


(2.4-81) 


VI = cosa, sa 2 

The actual sensitive axis of the equivalent vector seeker is 

determined by k = i x j whereas the assumed sensitive axis is 

coincident with I x 1, . The deviation between these axes is 
a b 

evaluted from 


cos a c 


(1 xl, ) . k 
a b 


'SA~|l a xl b | (2.4-82) 

From the projected views of Figure 2-24b 


it is simple to write 


1 = 
a 


■*"b a sa f 


— + o' sin n , 
l sa.^ 1 


sin r/ 2 i + 


- + a s a x cos „ l¥ 
- + a cos n 2 k 


(2. 4-83) 


The angles and p 2 are independent random variables, each 
uniformly distributed on the interval 0 to 2 tt. This recognizes that 
sensitive axis uncertainty for a single instrument is specified 
merely as an angular uncertainty from a nominal scribed line. 

If terms of second order smallness are 
retained when evaluating 2. 4-82 by use of 2. 4-83, the following 
results 


“SA 


Y 2 

{ a 

! sa i 


2 2 2 

cos m+a? cos r } 0 - a „„ a n sinrj 1 sin T7 0 > 
1 sa 2 2 sa^ sa 2 1 2J 

(2. 84) 


75 



In terms of mean square values, cr?. = a + a 

n SA saj^ sa 2 

which leads to the simple result for equal quality instruments in 
a null seeking arrangement. 


a 


2 

SA 


2a 


sa 


(2.4-85) 


Before proceeding to derive a g^ for other 
situations, let us examine an important distinction between the 
problems represented by Figures 2-22 and 2-24a. If a ^ is con- 
sidered on a deterministic basis, the value of is uniquely 
determined. In contrast, for a fixed value of a there is an in- 

S3- 

determinacy in the evaluation of Q'g^' This may be explained 
graphically as well as by writing the mathematical expressions 
for- agA . 


Consider i, j and k in Figure 2-24a as a 

fixed reference coordinate set. The 1 , 1, and (1 x 1, ) triad 

a b a b 

is considered as a rigid member free to pivot about its origin 
and, therefore, rotate with respect to the fixed triad. This 
motion is subject only to the constraint that the tips of 1 and 1^ 
remain in contact with the edges of the unit-height cones of half- 
angles a . If we assume the half-cone angles are identical, 

S3. __ 

each vector 1 and 1, may traverse a full conical locus while 
satisfying this constraint. One of the possible motions results 
in the line (1^ - 1 ) always remaining parallel to (j - i ) but 
varying in separation from zero to a . The angle ttg^ takes 
on corresponding values between zero and some maximum. The 
fact that dn. has a range of values for a fixed value of a is the 

bA S3. 

indeterminacy mentioned earlier. 

For an explicit solution of this determi- 
nistic problem, one must solve a three-dimensional geometry 
which mi ght, be described as a four-bar linkage mapped onto a 
sphere. The four links are arc-lengths subtending the angles 



a , tt/2, a and n/2 in sequence (See Figure 2-24c). The angle 
rj is the independent variable required to make the problem com- 
pletely deterministic. The reader is cautioned at this point that 
an explicit solution is afforded by the introduction of ri only be- 
cause the half cone angles a are assumed identical. 

& sa 


The mathematical expression for «g A is 


cosn SA = d a x V* (ixj) 

= <i a - i> < V j ) - U a - 1 > d b * i) 2.4-86 
Using 2. 4-81 with 

C os a SA = COS 2 a ga -• ^(I a .3)( 1 b - i) 2.4-87 

The unevaluated dot products of 2. 4-87 
are geometrically interpreted as the diagonals of the four-bar 
linkage mentioned above. These diagonals separate the four- 
bar linkage into the, two spherical triangles of Figure 2-24d. 
Application of the law of cosines yields the results 


1 .j = sin a cos r] 2.4-88 

ci sa 

1. . i = -sin a cos n 
b sa 

Substituting 2. 4-88 into 2. 4-87 and reducing for the case of 

small a „ 
sa 


a 


a 


SA T V “ “ sa Sin 11 2.4-89 

which has a maximum value of A = J'l a 

SA v sa 

There is another locus which satisfies 

the constraints 2. 4-81. In this case the lines (1, - 1 ) and 

— — b a 

(j - i) do not remain parallel but intersect. It happens that this 

results in the same maximum value for Therefore, for a 

nulling measurement with two sensors of deterministically equal 
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uncertainty., the equivalent sensitive axis uncertainty is given by 


0 y 2 a 

SA v sa 


2. 4-90 


The foregoing discussion relates to an 
admittedly restricted situation. Yet it is important to recog- 
nize the indeterminacies which may result from considering 
single-axis conical error loci. The concept of evaluating pa- 
rameters of an equivalent "vector seeker" suppresses these 
indeterminancies effectively. 

For a two-instrument magnitude measure- 
ment, an orthogonality constraint is required to describe the 
sensitive axis uncertainty as an equivalent rotation. In this case, 
the equivalent sensitive axis uncertainty is precisely the uncer- 
tainty in the best estimate of the individual sensors' axis loca- 
tions. For identical instruments a OA = a 

SA sa 

The equivalent "vector seeker" for a three- 
instrument magnitude measurement has its sensitive axis in the 
geometrical center of the triad. As in the earlier null-seeking 
situation, a triple constraint of the form 2. 4-81 does not yield a 
unique value of ag A even in a deterministic situation. Yet we 
know from Equation 2. 4-10 that these three constraints do fix a 
unique value for the magnitude of the whole angle rotation be- 
tween the two triads. The discrepancy involves the location of the 
axis about which the whole angle rotation develops. If the rota- 
tion happens to occur about the sensitive axis of the "vector 
seeker, " then <*g A = 0. However, if it develops about any axis 
normal to the sensitive axis, ag A takes on its maximum value. 

On a mean square basis, the equivalent sensitive axis uncertainty 
for a magnitude measurement with three sensors is given by 


a 


SA 


2 a sa 


2. 4-91 


The reader is reminded that the result 2. 4-91 is valid whether 
or not the sensitive axes are mutually orthogonal. We have, 
assumed, however, that they are. 



2. 4. 2. 5 Summary of System Measurement 

Capability 

At the beginning of Section 2. 4, the fol- 
lowing question was posed: "With what accuracy can the 

measurements made with sensors integral to a particular 
inertial system indicate the direction of a vector? " We review 
now the key results obtained in response to this question. 

The evaluation of the measurement capa- 
bility of a system depends upon computation of the statistics of 
the parameter a . Figure 2-12 shows the simple geometric 
interpretation of a as the angle between the indicated and actual 
vector directions. This single parameter serves for quantita- 
tive comparison of various systems regardless of the instruments 
or measurement technique under consideration. It has been shown 
(Equations 2.4-48 and 2. 4-49) how a is related to uncertainties of 
the instruments employed by a system and, further, that the 
values of a are a useful means for comparing the self-contained 
measurement portion of various alignment schemes. * 


The instrument uncertainties which affect 
alignment measurements fall into two classes, namely those 
related to knowledge of sensitive axis locations and those 

related to knowledge of transducer behavior (a,p). In addition to 
depending on the magnitude of individual instrument uncertainties, 
uncertainties for systems of sensors depend upon both the instru- 
ment configuration and the technique of measurement. The deter- 
mination of a from an d a^, however, is independent of the 
instrument configuration and measurement technique. 

The mean- square error contribution of 
2 

. sensitive axis uncertainty, <*g^, is presented in Table 2-2. En- 
tires in this table were analyzed in Section 2. 4. 2. 4. Because 
our definition of a does not admit negative values, a will 

S3. S3 

have an asymmetric probabili ty de nsity function with mean value 

2 

other than zero. This is why ^g^ is written in terms of mean- 


* The reader is reminded that a , therefore, relates to the 
maximum achievable accuracy, not the minimum. 



square values rather than variances. 

The mean-square error contribution of 
2 

transducer uncertainty, o^, is presented in Table 2-3. The 
transducer-limited ability of an inherently angle sensitive device 
to measure the direction angle of a vector is obtained directly 
from its specifications, a ^ denotes this uncertainty. It is mea- 
sured in a plane containing the vector and the instrument's sen- 
sitive axis. The complete specification of a vector's direction 
only requires the measurement of two direction angles. The 
corresponding measurement planes intersect in the vector and 
it is assumed they are orthogonal. The angular uncertainty, 
a within which two such sensors indicate the vector's 

/ ~2 2 ^' 2 

a , + a t . (In Table 2-3 a . _ 

*1 2 t l 

~2 2 

af- = af. . ) The entries in Table 2-3 pertinent to sensors 
whose output is proportional to the magnitude of the vector mea- 
sured, summarize the efforts of Section 2. 4. 2. 2 and 2. 4. 2. 3 in 
converting magnitude uncertainties into equivalent angular un- 
certainties. 

Tables 2-2 and 2-3 are meant primarily 

to help evaluate the mean- square direction -indicating error of 

a system. Knowing the measurement technique to be employed, 

2 

th e va lue of obtained with the help of Table 2-2 and the value 

of obtained with the help of Tab le 2 -3 are combined by the sim- 
1 2 
pie sum of Equation 2. 4-48 to yield a . This sort of procedure is 

the essence of the "vector seeker" concept. All measurement 

techniques and sensor descriptions are reduced to an ability to 

seek the direction of a vector. An evaluation of other statistics 

of the vector seeker, as for example the mean value of a (See 

Equation 2.4-48), requires further reduction of the equations 

for a derived in the text. 



As a somewhat trivial example of the use 
of the tables, we calculate the maximum accuracy to be expected 
in finding the direction of g with particular accelerometers. From 
a standard catalog, we chose at random an accelerometer with 
specifications (as interpreted by this author) 

a =0.1 mr 

sa 

a = 5x10 ^g 

3 . 


^aj = 1x10 4 gjg 


By , simple algebraic manipulation we find 


a = 0. 158 mr for null- seeking 

rms & 


a 


rms 


= < 


0. 163 mr max 


0. 141 mr min 


for a triad 


In addition to providing a calculational 
aid. Tables 2-2 and 2-3 yield certain fundamental comparisons 
between various measurement methods. The results of Table 
2-3 have been largely discussed in Section 2. 4. 2. 2. One remaining 
comparison is that for measurements made with angle- sensitive 
instruments to be equivalent to those made with magnitude -sensitive 
instruments, then either 
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depending upon whether null-seeking or direction-seeking are 
employed respectively. 


Equation 2. 4-48 shows that on a mean- 
square basis the contribution of sensitive axis uncertainty to 
a system's vector indicating capability is every bit as important 
as the contribution of transducer uncertainty. Reduction of the 


magnitude of is, thereby, an i mpor tant objective. Table 2. 3 
shows that a reduction of 25% in a jf may be had by using two 


instr uments rather than t hree . An alternative means of reducing 

2 2 

or 0 . is by minimization of a 

bn S3. 


The ultimate goal of alignment is to orient the 
guidance sensors (gyros and accelerometers). Yet frequently 
this is achieved indirectly by aligning a set of fiducial axes scribed 
in the IMU's stable member. In any situation where a set of 
scribe-marks define the prime reference frame of a system, the 
relative orientation between the sensitive axes of the inertial sen- 
sors and reference coordinates must be established by a carefully 
controlled instrument mounting procedure and/or a calibration pro- 
cedure. Knowledge of this relative orientation is necessary in order 
to transform data from the measurement frame defined by the sensi- 
tive axes to the scribed reference frame. It is precisely such a 
data transformation requirement which introduces sensitive axis 
uncertainty problems. Let us pursue this by assuming the role of 
an observer standing on the platform. If an inertial vector is mea- 
sured for alignment, the inertial sensors can determine their orien- 
tation relative to the vector (e. g. , the accelerometers locating g). 
But unless we know our orientation (i. e. the scribe-mark locations) 


relative to the instruments, this information is useless in establishing 
our orientation relative to the vector. describes the uncertainty 

in knowledge of our orientation with respect to the instruments. But 
of what importance is our orientation (scribe-marks) to the guidance 


function? In a word, none. When measurements of the inertial 


sensors are used for alignment as well as for guidance, the system's 
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prime reference frame should be affixed to the sensitive axes 
of a set of inertial sensors. The location of the reference 
frame with respect to the structure is no longer of interest, 
except to know that it is described by a fixed rotation. In just 
a very few words, by defining our reference coordinate frame 
coincident with a set of instrument axes we set the sensitive 
axis uncertainty for those instruments to zero. 

Lest anyone promote the setting of a 
reference frame according to an instrument frame as a cure- 
all for sensitive axis uncertainty- problems, there are many 
limitations. First, it requires a triad of sensors to do a job 
of which a pair of sensors may be capable. Second, unless one 
triad of instruments can be used for all alignment measurements 
(at least two required), the problem of relating two sets of sensors 
to one another is very similar to relating each set to an arbitrary 
reference frame. Third, external information supplied to the 
system must be coordinatized via the platform axes seen by the 
external sensor. This requires the sensitive axes to be cali- 
brated in platform reference coordinates once again. More 
particularly, when non-inertial inputs are used for alignment, 
the use of scribe-mark coordinates are necessary to relate the 
sensitive axes of inertial and non-inertial sensors. 

In view of these drawbacks, it is best stated 
here only that sensitive axis uncertainty 1) is a significant source 
of error, 2) can be reduced by careful mounting and calibration 
procedures, and 3) can be eliminated for one set of sensors by 
the proper choice of reference coordinates. 

2. 4. 3 Composite Misalignm ent For Two Vector Me asure- 
ment 

Thus far the ability of an IMU to independently indi- 
cate the direction of, each of two non-collinear vectors has been 
established. This ability is expressed for each vector by the un- 
certainty angle a. It remains to develop the- magnitude of the total 
IMU misalignment resulting from imperfect measurement of the 
vectors. 
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The directions of the two vectors (V^, Vg) being 
measured and the loci of possible independently indicated di- 
rections are shown in Figure 2-25. 9 is the angle between the 

measured vectors, corresponding for example to the co-latitude 
angle in a gyrocompassing system. The dependent reference 
frame being aligned is essentially fixed to the unit vectors (ly )^ n( j 

and (fy )^ nc j- It may be represented by the non -orthogonal 
2 


coordinate 


set composed of (ly^ ind> (ly^^, 


and the unit nor- 


mal to them, (i" n ) ind - The frame being aligned to is similarly 

fixed to the coordinate set \, T , F Tr , T . 

V 1 V 2 n 

It is shown at the end of Section 2. 4. 1 that, provided 


(1 V 1 ) ind x ^V^ind 


( 1 v x x 1 v 2 ) 


sin 9, 


the evaluation of the magnitude of the whole angle rotation, $, for 

the non- orthogonal coordinate set is the same as the evaluation 

for an orthogonal set. The provision that the angle between 

(T. r ). , and (T tt ) . , be 9 is more a naturai than an artificial 

V, md V, ,md 

constraint. As noted recently, 9 is the co-latitude angle in 
gyrocompassing system. . In master-slave configurations, 9 
is effectively measured by the master and supplied to the slave 
for control purposes (9 is implicitly tt / 2 for vertical-azimuth 
situations). Lack of any constraint between the allowed indica- 
tions merely renders the concept of a total misalignment mea- 
ningless. In that case the realizable platform misalignment 
angles assume the infinity of values related to all combinations 
of points on the indicated loci. 

The angles a ^ and a 2 will be recognized as equiva- 
lent to the (f>^ of Equation 2.4-1 so that $ is obtained from 
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Figure 2-25 Geometry of Two-Vector Composite Misalignment 
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* 1 / 2 , 2,2 
~JTJ 1 " 2 “ 3 

where a g is defined by Figure -2-25. To calculate a. 


2. 4-10 


COS O!^ = (1 ) * (1 ). , 

3 n n md 


e r Vl xl V 2 ] [ u V 1 'ind x ' 1 V 2 'ind] 
sin 0 ' sin 0 

= - - 9 - < COS a cos a 9 - [f . (1 ). ] 

sin 0 L 1 2 V 1 V 2 md 


(it;- x (i,, ). 


*- 1 V 2 * (1 V 1 ) ind^ 


from which 


— ^9— | a ? + a\ ~ 2 cos^ 0+2 

sin 0 ] 1 2 


[>r 


(1 V 2 ) in^l 1 V 2 - ^V^indj 


2. 4-93 


Except in the vicinity of 0 = 0° or 0 = 90° the approximation 


ly i ’ <lv 2 )ind “ 1 V 2 * (lv i )ind “ C ° S 9 


2. 4-94 


is reasonable. Cases of 0 near zero are uninteresting but cases 
of 0 near 90° are important. Figure 2-24 and its associated 
mathematics substantiate that the value of resulting from use 
of the approximation 2. 4-94 at 0 = 90° is the correct upper bound 
on and therefore also on$>. (The value of a g is only bounded 
for fixed a a 2 and 0; it is not unique. ) 

Introducing Equation 2. 4-94 into 2. 4-93 and intro- 
ducing that result into Equation 2. 4-10, the following is obtained: 


= 2 + a 2 ^ + csc ^ 


2. 4-95 
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2 2 

-The mean-square angular uncertainties, and are n °t g e ' 

nerally equal. For cases where they may be related by a simple 

— 2 ~ 

num er ical factor n, ascribe the literal a to the smaller of 
2 2 

a i and such that in the expression 2. 4-96 

* 2 = [-|-(n 2 + 1)(1+ esc 2 0)] a 2 2.4-96 

the quality factor n is always greater than unity. The square- 
root of this equation is plotted in Figure 2-26 as a function of 
n and 0. 

A word regarding the "quality factor" n is in order. 

It is admittedly a shaky proposition to relate complex systems 
by simple numerics, but only from the viewpoint of mathemati- 
cal rigor. From the viewpoint of practical utility, the use of 
simple numerics for comparative purposes is commonplace. 

That we have chosen the ratio of root-mean- square direction- 
indication errors as a figure-of-merit is a natural result of the 
preceeding developments. 

It may be seen from Figure 2-26 that at 0 = 90° the 
total rms misalignment is approximately 98 per cent attributable 
to the poorer measurement when n = 5. In fact, when one mea- 
surement is only twice as good as the other at 0 = 90°, the poorer 
measurement accounts for almost 90 per cent of the total error. 
Only when 0 is significantly less than 90° does the total error 
reflect a meaningful component from the more accurate measure- 
ment. 

2.4.4 Information Transferred From An External Sensor 
or System 

Section 2. 4. 2 was concerned with the ability of a 
set of sensors to indicate the direction of a vector input. One 
of the characteristics of moving base phenomena is that frequently 
the measured vector differs in direction and/or magni- 
tude from the desired vector. If disturbances causing changes 
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in the desired vector's direction are of sufficient intensity, com- 
pensatory information may be supplied in an attempt to reduce 
the discrepancy. Considering that compensation procedures are 
rarely perfect, a residual angular uncertainty ( 0 ) between the 
direction of the desired and compensated actual vectors (See 
Figure 2-27) must be accounted for. The resulting total uncer- 
tainty in finding the direction of a desired vector input is assigned 
as Y. 

From Figure 2-25, 

Y* = or 2 + 0 2 -2 a 0 cos rj 2. 4-97 

where r\ is uniformly distributed. Using the results of Appendix 
B, Equation B-13 and B-16, 

T = — (cH-0 ) t (k 2 ) • k2 = — £ 2, 4-98 

"■ (a+pr 

f = a 2 + 0 2 2.4-99 

a was the subject of Section 2. 4. 2. The evaluation of 0 , since 
it is strongly related to specific schemes applicable to the 
moving base alignment problem, is deferred to Chapter III. If 
0 is non-zero when the results of Section 2. 4. 3 are used, a 
simple substitution of yfor a is rigorously correct. 

2. 5 Summary 

This chapter has attempted to develop a unified theory of 
alignment whose purpose is to lend coherence to the analysis of 
numerous techniques proposed as solutions of the moving base 
alignment problem. It is founded on the two considerations common 
to all alignment procedures: coordinate frames and measurements. 
This serves to adequately bound expected alignment accuracies. 

All alignment procedures are viewed ks the transfer of orientation 
between two coordinate frames by the measurement of at least two 
non-collinear vectors. Separate techniques vary only in the me- 
chanization of particular measurements. 
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Figure 2-27 Geometry Relating <yto a and /3 
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Coordinate frame rotations are handled as whole-angle 

/ 

rotations, the latter being developed in terms of all the more 
usual rotational parameters. Certain interesting results are 
derived concerning the magnitude of the whole-angle rotation 
in the case of Euler angles and principal direction cosines. 

The latter angles are also shown to be an especially convenient 
representation for problems involving non- orthogonal coordi- 
nate sets. 

Measurements are discussed from the viewpoint of a 
system of single-degree-of-freedom sensors attempting to 
locate the direction of a vector input. Both magnitude- sensi- 
tive and angle sensitive instruments are considered. In the 
former case, bias and scale factor uncertainties are the assumed 
errors. In the latter case, the uncertainties are inherently 
angular. Sensitive axis location uncertainties are included for 
both cases. The measurement error of any system is reduced 
to an angular uncertainty in indicating the direction of a vector 
by converting the uncertainties of magnitude sensitive instru- 
ments to their angular equivalents. This conversion is carried 
out for all practical combinations of sensor-to-vector orienta- 
tion. An important result of this conversion is the establishment 
of definite guidelines to minimize magnitude and direction indica- 
tion errors as a function of instrument uncertainties and the 
orientation of the vector with respect to the sensors. 

Deviation of a vector input from a desired nominal direc- 
tion due to base motion is expressed in angular terms and inclu- 
ded with the basic measurement error. The magnitude of the 
total platform misalignment resulting from independent mea- 
surements of two non-collinar vectors is calculated as a function 
of the angular uncertainties about the nominal positions of each 
of the vectors. 

Application of the unified analysis to specific alignment 
techniques is accomplished in Chapter III. 





CHAPTER III 


ANALYSIS OF ALIGNMENT TECHNIQUES 


3. 1 Introduction 

The preceding chapter considers all approaches to the 
solution of the moving base alignment problem as a class of non- 
dimensional vector direction indication problems. In contrast to 
this, the majority of the literature concerned with aligning inertial 
systems on a moving base is an exposition of distinct "techniques. " 
The abstractness of the former approach is reconciled with the lat- 
ter approach in this chapter. This is accomplished by individually 
discussing several techniques within the analytical framework of 
Chapter II. 

The most significant techniques appearing in the literature 
(and those considered in this chapter) are: 

1) Vertical Indication 

2) Gyrocompassing 

3) Star Tracking 

4) Fix Monitored Azimuth 

5) Vector Matching 

6) Gimbal Angle Matching 

7) Optical Comparison 

Each of the foregoing techniques makes measurements on a dif- 
ferent vector quantity. This is perhaps the clearest distinction 
among them. Once the vector associated with a given technique is 
identified and the measurement sensors are defined, the limiting 
ability of the sensors to indicate the direction of the vector is ob- 
tained from the results of Chapter II (as a quantitative evaluation 
of the angle a). We also establish the significance of the vector 
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relative to a coordinate frame of interest and express the degree 
to which it maintains its significance in the face of perturbing en- 
vironmental influences by evaluating the angle ft. The total vector 
direction indication error, y, follows directly from Equation 2.4-97 

The techniques by which alignment is effected fall into two 
distinct categories: Direct Measurement Techniques and Transfer 
Techniques. Recalling the definition of alignment, angular coin- 
cidence is sought between the dependent reference frame (Section 
2. 2. 2) and the independent reference frame (Section 2. 2. 1) of Figure 
2-1. This is accomplished by measurements on at* least two non- 
collinear vectors. In direct measurement techniques, equipment 
associated with the dependent reference frame measures vector 
quantities which have a unique geometrical relationship to the in- 
dependent reference frame. This geometrical relationship con- 
stitutes a priori information which must be available to the depend- 
ent reference frame. A simple example is the measurement of g 
coupled with the knowledge that it lies along one axis of a local 
geographic coordinate frame. 

In transfer techniques, equipment associated with the de- 
pendent reference frame measures vector quantities which do not 
necessarily have an a priori geometrical relationship to the in- 
dependent reference frame. Rather, an intermediate (Section 2. 2. 3) 
or transfer (Section 2. 2. 4) reference frame simultaneously 
measures the same vector quantity as the dependent reference 
frame. This second measurement serves to establish the geomet- 
rical significance of the vector with respect to a mechanized 
representation of the independent reference frame. Once deter- 
mined, of course, this information must be communicated to the 
dependent system. The vector, as a common sensory input to 
both systems, is serving only as an information transfer medium. 

Of the measurement techniques listed earlier, the first four are 
in the direct measurement category; the remainder are in the 
transfer category. 

Prior to examining individual techniques, we call attention 
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to one further distinction. Indicating the particular direction known 
as the vertical is well enough solved in both theory and application 
as to be the most widely acceptable measurement for moving base 
alignment. In fact, because of the tendency to rely upon the ver- 
tical to provide one of the two required vector inputs, the moving 
base literature is largely oriented to the single vector problem of 
azimuth determination. For this reason, many of the techniques 
discussed later in this chapter are intended only for azimuth align- 
ment of a system which is already correctly oriented with respect 
to the vertical. 

3. 2 Vertical Indication 

We consider here the fundamental limitations in locating 
the vertical direction by means of specific force measurements. 

The specific force sensor employed for vertical indication on a 
moving base is the accelerometer. Devices such as a spirit level 
are often employed in a more benign environment but do not pos- 
sess the combined dynamic range and accuracy of which acceler- 
ometers are capable. We are interested in quantifying how well 
the accelerometer can measure a specific force vector and how 
well this vector 1 s direction represents the vertical direction. 

The independent reference frame for vertical indication is 
an earth fixed reference frame. This reflects the unique geomet- 
rical relationship between the vertical direction at a point and the 

13 

earth. According to Wrigley 

"The Vertical at a point on the Earth' s surface is 
defined as the local direction of the force of gravity 
as indicated by a plumb bob hanging from a base 
that is stationary relative to the Earth. The hori- 
zontal plane is defined as the plane perpendicular 
to the vertical, or as the surface of a free liquid 
under the influence of gravity alone. " 

Note that although the vertical is intimately related to 

a physical vector, it is not itself a vector. It is merely a signed 

direction. This is as stated in Section 2. 3. 1: when a vector is the 

input to the alignment problem, only its direction is of ultimate 

interest. There is no requirement to explicitly compute the 


97 


magnitude of the gravity vector and, therefore, no requirement to 
perform a complete vector measurement. 

In considering vertical indication from a carrier operating 
above or below the earth's surface, the local vertical is under- 
stood to imply the direction of the extended plumb bob vertical 
passing through the center of measurement of an inertial system. 
Additionally, we assume that deviations in the direction and mag- 
nitude of gravity (g) are generally negligible within the volume of 
space occupied by a carrying vehicle. For an analytic gravity field 
non- parallelism is only on the order of 0. 01 arc second per foot 
of separation in the horizontal plane. By comparison, gravity 
anomalies may cause parallel gravity vectors to appear at points 

/ 17 

on the earth 1 s surface within a subtended angle of 15 arc seconds. 

Of far greate/ significance to the vertical indication prob- 
lem is the difference between the "dynamic" (or apparent) vertical 
and the "static" (or true) vertical. The true and apparent verticals 
are defined respectively by the directions of gravity (g) and total 
specific force (f). Depending upon one's point of view, either the 
apparent vertical is an imperfect representation of the true ver- 
tical or the true vertical is a special case of the apparent vertical 
occurring when the measuring systern is at rest with respect to the 
earth. Since only the true vertical has a unique geometrical mean- 
ing with respect to the earth, we adopt the former point of view. 

That is, by severing the connection between the system and the 
earth, we invite base motion to destroy the uniqueness of the quan- 
tity we are able to measure. Error- free accelerometers are only 
capable of indicating the direction off; they require additional in- 
telligence to deduce the direction of g from that of f. In essence, 
we must now process the moving base measurements in order to 
obtain an equivalent fixed-base indication. 

Relative to the format of Chapter II, the angular discrep- 
ancy between g and f is geometrically defined by |3. (See Figure 2^-27 
reproduced here as Figure 3-1, but with the vector designations 
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Geometric Quantities Relevant 
to Vertical Indication 
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appropriate for vertical indication. ) The problem of deducing 
the direction of g from f is equivalent to forcing /3 to zero. The 
limited instrument capability to accurately indicate the direction of f 
is given by the appropriate angle a (relative to the instrument con- 
figuration employed) 

3.2.1 Evaluation of ft for Vertical Indication 

A small angular difference (/3) between the true and 
apparent verticals is described by 

g . Ig * £l 

l?l Ifl 

For most conditions of interest. 



and _ _ 

o -IjjJill 

p ” g 

3. 2-1 

from which it is seen that the horizontal component of f deflects 
the vertical at the rate of 1 minute of arc for each 290 micro- g. 

A mathematical relationship between g and f is 
obtained in convenient form by expressing the inertially referred 
acceleration of a point as the vector sum of component acceler- 
ations referred to several coordinate systems. Using the four co- 
ordinate frames and position vectors defined by Figure 3-2, f at point 
D may be written as 

7 - G- Pj 2 r id 

f “ G - [p j Rjg + P E u IE x R ed 1 

"|pE R EC + P E W EC X R CD + 2u EC x p C R CD + U EC x(u EC x R Cd|] 

"[ 2u IE x p E R ED + W IE X ^ W IE X “ P C R CD 3 ’ 2-2 

* See Appendix D, Reference 13 for a detailed mathematical 
development. 
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I is center of an inertial system - no acceleration relative 
to fixed space 

E is center of the earth; origin of axes fixed in the earth 

C is center of gravity of the carrying vehicle; origin of 
axes fixed in the vehicle 

D is center of dependent reference frame; origin of axes 
fixed in the inertial measuring unit 

**ID = ^IE + ^EC + ^CD 


Figure 3-2 Vector Representation by Components of the 
Position of the Center of Measurement of an Inertial 
Measuring Unit in Motion with Respect to the Earth 
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where 

P i is the time derivative operator in the subscripted 
coordinate frame 

W EC * s t ^ ie an S u ^ ar velocity of the vehicle about its 
center of gravity measured with respect to the 
earth 

w JE is the angular velocity of the earth with respect to 
an inertial frame 

G is the gravitation at point D 

Equation 3.2-2 is equally valid for vehicles opera- 
ting below, on or above the surface of the earth. Figure 3-2 
merely depicts an extra-terrestrial vehicle for convenience. 

In assuming that f is the specific force measured 
by the IMU's accelerometers, we ignore a slight refinement oc- 
casioned by the required pendulosity of the instruments. A more 
complete characterization only results in additional terms of 
negligible magnitude. 

As an initial simplification of Equation 3. 2 -2, 11. c terms 
included in the first bracket make no significant contribution to 
f. They will, therefore, be dropped. Furthermore, the definition 
of the observed local acceleration of gravity, "g, at point D 

g * G - x (Wj E x Rjjjj)) 3.2-3 

allows two terms to be combined. These simplifications result in 
the following expression for f 

f = g " [p e R E c + Pe w EC x R CD + 2u EC X P C r cd 

+ w EC x (w EC x R cd ) + P c R cd + 2w IE x p E R ed ]] 3 2-4 

Bracketed in Equation 3.2-4 a;.e the terms which corrupt 
the physical observable representing the vertical. We are inter- 
ested in their horizontal components. The physical sources of the 
interference are motions of the vehicle relative to the earth and 
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dynamic behavior of the vehicle as a non-rigid base. Notice that 

static alignment accuracies between the C and D frames are not 

important. In the absence of system motion with respect to the 

earth, f = g, whereas, in the face of significant base motion, the 

magnitude of j8 (Equation 3.2-1) is sufficient to require compensatory 

modification of the accelerometer's raw outputs. Certain terms 

of Equation 3. 2-4 are known, in fact, as "acceleration compen- 
17 

sation terms. " When compensation is employed. Figure 3-1 
is still accurate for error determination provided /3 now 
represents the residual uncertainties for the compensation 
process. 

The compensation requirements for individual appli- 
cations are derived in a straightforward manner from Equations 
3.2-1 and 3. 2-4. However, the requirements are too strongly 
related to parameters such as desired accuracy, vehicle config- 
uration, vehicle speed and expected maneuvers to permit general 
quantitative conclusions here. We can discuss, however, the al- 
ternate means of compensation. 

Compensation of the output of a sensor (or some 
deterministically derived function thereof) is basically of two 
types: 1) signal modification by filtering; and 2) signal modification 
by the addition of an externally derived signal. Depending on the 
character of the base motion, these techniques may be employed 
separately or simultaneously. The attraction of signal filtering is 
that it may be mechanized on a self-contained basis without addi- 
tional sensory requirements. Unfortunately, it is limited in its 
ability to discriminate between signal components on any basis 
other than frequency. Externally derived compensation signals, 
while imposing the requirement for additional sensory equipment, 
are capable of accurately discriminating between like frequency 
signals on an amplitude basis. Filtering also may be applied directly 
to the output of external sensors in order to reduce the noise 
associated with these measurements. 

The frequency at which g is observed depends only on 
the relative motion between the dependent reference frame and 
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g. By Equation 3. 2-4, angular motion of the clop-.mJont frame does not 
affect f. Therefore, the frequency modulation of g by platform 
motion is reflected in the same fashion by all terms of f. This 
permits us to arbitrarily choose to examine base motion inter- 
ference from a frame in which "g is a zero frequency quantity. 

(Note that this restricts only two degrees of freedom of the plat- 
form. ) 

As to motion of the carrying vehicle relative to the 
earth. Brock* states that the filtering problem is greatly aided if 
the vehicle is controlled so that the velocity (vector) between ma- 
neuvers fluctuates only slightly from some nominal value. An 
examination of Equation 3, 2-4 shows tnai this is indeed a help in re- 
ducing bias accelerations. Yet for the very condition of constant 
velocity, a bias Coriolis acceleration (2 Uj E x pgRgp) deflects 

_ 9 

the vertical at the rate of 2. 29 x 10 sin X minutes of arc per 
mile per hour of horizontal speed (X = latitude). For a ship at 
high latitudes this might approach a value as high as one arc min- 
ute. For a subsonic airplane it might easily approach 10 arc 
minutes. This term will generally have to be compensated, re- 
quiring thereby sensors capable of providing a measure of the 
vehicle velocity p^R^^. 

A measure of the high bias accelerations generated 

2 — 

by vehicle maneuvers is given by the term p E R E £ eva l ua t e d at 
constant speed. The radial accelerations at constant speed are 
plotted in Figure 3-3. Note that a man walking at 5 miles per hour 
around a circular path of one mile radius experiences a constant 
vertical deflection of one arc minute. While these high maneuver 
accelerations make vertical indication difficult, they are indis- 
pensable to the vector matching technique of Section 3. 5. 

An example of an easily filterable horizontal accel- 
eration is that due to vibrational modes of the vehicle. Given a 

2 — 

suitable period of observation, physical constraints prohibit PqRqq 
from making any net contribution to f. This is true even in a 
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violent environment such as a missile under the wing of an air- 
plane. The time required to discriminate against this (and any 
filterable quantity) is determined by the vibrational spectral 
characteristics, the nature of the filter, and the amount of noise 
suppression required. 

Schuler tuning is, of course, a unique form of fil- 
tering whereby the vehicle's tangential acceleration becomes an 
aid in tracking the vertical rather than a quantity to be discrim- 
inated against. This concept and the attendant damping problems 
are thoroughly considered in the literature. 

3. 2. 2 Evaluation of a for Vertical Indication 


Measurement of the specific force vector's direction 
using single- degree-of- freedom accelerometers requires either 
two or three instruments. With two accelerometers, they must 
be servo- controlled to enable orientation of their sensitive axes 
in the null plane. This represents by far the most common mech- 
anization of vertical indication. From Tables 2-2 and 2-3 (with 


so that 


4 - ^ f* o/ g ) 2 

?-[*£♦ fS a ] 


by Equation 2. 4-48. 


For three accelerometers, it is beneficial but not 
required to specially orient the sensitive axes. The tables yield 


2 3 2 

= — a 
2 sa 





cos^i] 


3. 2-6 
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In Equations 3.. 2- 5 and 3.2-6, a is the mean- square 

sa 

sensitive axis uncertainty of each instrument expressed in radian 
measure to the second power, C a ( - ) is the standard deviation of 
accelerometer bias expressed in equivalent "g's, " and cr^ a is the 

standard deviation of accelerometer scale factor (a dimensionless 
quantity sometimes expressed in units of gig). Notice that with 
either two or three accelrometers, vertical indication accuracy 
is inherently bias limited. 


3. 3 Gyrocompassing 

I 1 7 

The interpretation of the term "gyrocompass" (the pre- 
cursor of gyrocompassing) has remained unambiguous, that being 
a pendulous gyroscopic element arranged to track the horizontal 
component of the earth's daily rotation The term "gyro- 

compassing, " on the other hand, has become subject to much 
looser interpretation, to the point that it more connotes the gen- 
eral process of aligning by means of an angular velocity vector 
than it does the specific process of aligning to by a particular 
mechanization. We restrict our attention here to the vector 
and investigate the fundamental limitations in establishing its 
direction with gyroscopic sensors mounted on a moving base. 


The unique geometrical relationship between the direction 
of Ujg and the earth's polar axis dictates that the independent ref- 
erence frame be an earth- fixed reference frame. This is the same 
frame employed for vertical indication. Base motion relative to 
this frame, it will be shown, generally can cause severe bias de- 
flections of the measured angular velocity relative to Uj^,. 

Prior to evaluating a and 13, let us briefly consider a prac- 
tical example of the importance of determining system alignment 
accuracy relative to instrument performance. The majority of 
operational situations requiring alignment of a dependent frame in 
a mobile environment will employ the dependent frame, once aligned, 
for a relatively short period of time. In such situations, gyros with 
only moderate drift performance may be capable of maintaining the 
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dependent reference frame with sufficient accuracy. The prac- 
tical advantages inherent in the use of such gyros are numerous. 

It is desirable, in such situations, that the alignment scheme not 
require a level of gyro performance better than that required for 
guidance. Thus, although gyros of sufficient performance to per- 
mit gyroscopic alignment may be within the current technological 
capability, such alignment methods must often be ruled out. The 
uniform presentation of performance parameters developed in this 
thesis is an aid in making such a decision. 

3.3.1 Evaluation of 6 for Gyrocompassing 

Referring to Figure 3 2 the angular velocity relative 
to inertial space sensed by the dependent frame gyros, * s 

expressed as the simple vector sum 

W ID * W IE + W EC + W CD 3.3-1 

An angular difference (/3) between the direction of and the earth's 
polar axis is described by 


sin jS 


l W IE x W ID 


U) 


IE 


w. 


'ID 


3. 3-2 


We must defer temporarily an explicit solution of 3. 3-2 for j3. As 
mentioned earlier, the possible base motion influence here is much 
more severe than that found in vertical indication. Gravity is a 
field force; its influence is felt even when removed from stable con- 
tact with non- force motion- parameters such as 1 ^ . One of the 
better advertised "problems" of angular rate sensing is that a 
westward speed of approximately 1040 cos X mph (X « lattitude) 
cancels completely (Equation 3 . 3-2 is then indeterminate. ) 
Rather than fearing the lack of a quantity to sense, one should be 
more concerned that a reduced magnitude vector input always 
increases the indication error (See Table 2-3), Hovorka* 5 concludes 
from the masking of earth rate by vehicle speed that gyrocompassing 
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is a slow-vehicle process. Independent of the rationale, this is 
a widely held viewpoint . 


To the extent that the dependent frame is caged to the 
vehicle, represents structural angular vibrations of the 

vehicle. If the structure cannot assume quasi- static deflections, 
<^CD must have zero-average value; if the structure does assume 
quasi-static deflections, the average value of w^^tends to zero 
over an extended period of observation. It appears, therefore, 


that can be successfully filtered from the gyro's output. 

This we will assume, with the understanding that it may only 
be a first approximation. Even if has a zero average value, 

the gyro output may well have a non- zero average indication of 
due to "kinematic rectification" ^ or "coning. "22,23 
The magnitude of a problem such as this is strongly related to 
design details and its examination is beyond the scope of this 
thesis , 


If we set to zero in 3.3-1 and assume that the carry- 

ing vehicle ig restricted to motion in the local horizontal plane, 
Eq 3.3-2 may be evalueated with the result 



where V is the vehicle groundspeed, is the northerly component 
of groundspeed and R is the distance of the vehicle from the center 
of the earth. Notice that even at constant speed and constant head- 
ing, ft may have a bias value which is non-filterable. 

The problem of forcing /3 to zero is equivalent to forcing 
V^j- to zero. This may be accomplished obviously by heading due 
east or west, or, in fact, by stopping the vehicle. It may also 
be accomplished by an external source of velocity .inf ormati on 
(such as Doppler radar, airspeed indicator or pitometer log) which 
provides a signal to compensate the gyro output.- Any one of these 
situations serves to place the measured vector collinear with its 
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meaningful orientation. The compensation process effectively 
re -couples the gyros to the earth. 

It is recognized that errors will exist in any auxiliary 
reference and that their nature will vary. Let us assume that the 
reference is subject to bias errors (a^) and sensitivity errors 
(Aa^). Bias errors may exist simply because speed is measured 
relative to a medium (air, water) which itself is in motion 
relative to the earth. Sensitivity errors will be sensor related 
but may correlate with the environment, as in the gain shift of 
Doppler radar with varying terrain reflectivity. 

When such a compensating reference is used, Equation 3.3-4 
becomes 


0 = 


(a Q + Aa 1 V N ) 

R(o;2 e + (a Q + A ai V) 2 /R 2 ) l/2 


3,3-4 


which, if the reference error coefficients have zero mean value 
and we neglect the second order effect of an error in the total 
velocity measurement, has a mean square value 


a 


p 2 = 


a. 


2 

+ CT A a] 

R 2 to ^ 
n IE 


V. 


N 


3.3-5 


It is noted that in developing this expression we have neglected an 
additional second order correction due to imperfect resolution of 
total measured velocity into north and east components . Details 
such as this must be considered at the design stage. 

In Equation 3. 3-5, °ag is the standard deviation of velocity bias 
expressed in the same units as (feet per second, miles per 
hour, etc. ) and , is the standard deviation of velocity scale 

factor (a dimensionless quantity). As an indication of the order of 
magnitude of these standard deviations. Reference 24 states that 
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the accuracy of modern Doppler radar approaches 0. 1%; o^ a is 

thereby implied of the order of 0. 001. Uncompensated northerly 
groundspeed deflects co ^ at the rate of 1 arc-minute per 0.3 miles 
per hour. In a 300 mph plane using 0. 1% speed measurement, 
fi has at least a one arc-minute value for northerly flight. 

3.3.2 Evaluation of a for Gyrocompassing 

Measurement of the direction of the angular velocity 
vector, requires two or three single-degree-of -freedom 

gyros. With two gyros, they must be servo -controlled to enable 
orientation of their sensitive axes in the null plane. From 
Tables 2-2 and 2-3 (with | w jq | » \ <* m \ ) 


a 


2 

SA 


a 


2 

T 



so that 


a 


2 £ a 


sa 




3.3-6 


With three gyros, the tables yield 

/cr 


2 3 2 „/ a 

a = 2 ff sa + 2 


°/“l E ) + ° Aa 


1 - 2 cos ip. 
1 \ i» 1 1 


3.3-7 


One of the more common system mechanizations places one gyro 
east, one north, and one either parallel or anti-parallel to the 
gravity vector. In this case, 3.3-7 becomes 


a 


_3_ 2 ^ 

2 ff sa + 



(1 - 


4 . 

cos X - 


. 4 . . 
sin X ) 


3.3-8 
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where A = lattitude. 


In these equations, a is the mean square sensitive axis 

S cl 

uncertainty of each instrument expressed in radian measure to 
the second power. 


a 


is the standard deviation of bias gyro 


drift expressed in the same units as w jg, and u^ a is the 

standard deviation of gyro scale factor ( a dimensionless quantity). 
Notice that with either two or three gyros, the directional indica- 
tion is bias limited. 


The total error, 7 , in indicating the direction of w j-g is 
described geometrically by Figure 3-4 (See also Figures 3- 1 and 
2-27) and the mean square value is given analytically as 


2 


7 



+ 




2.4-9 



Figure 3-4 Geometric Quantities Relevant to Gyrocompassing 


112 



3.3.3 A Special One -Gyro Situation 


12 17 

An important gyrocompassing system mechanization 2 * * 5 * 
relies upon accelerometer measurements to place the sensitive 
axis of a single gyroscope in the horizontal plane. The gyro axis 
is also adjusted in the plane until it detects no component of earth 
rate's horizontal component. By this process, the gyro's input 
axis indicates an aximuth line precisely 90° from true north. 

In the context of Chapter II, this is a single instrument 
nulling on a vector component of magnitude ( w jg cos X ). From 
Table 2-3, the transducer error is described by 


2 

Of rjn 



3.3-9 


This, however, is not the complete story. Whereas a sensor is 
insensitive to rotations about a vector, it is certainly sensitive 
to rotations about a component of a vector. A vertical error of 
0y n about a north- south axis causes the gyro to sense a compon- 
ent of earth rate's vertical component of magnitude (0y n Wj sin X). 
Equation 3.3-9 is corrected to read 


2 

Of rjn 



COS 



2 2 * 
0y n tan X 


3.3-10 


The total mean- square error in azimuth indication is given by 

r >2 2 2 ,.2 q o ii 

0 a = ff sa + q T + 13 3.3-11 

2 

wh ere a is the gyro's m ean -square input axis uncertainty, 

2 S3. 2 

a .p is given by 3 . 3-10, and /3 is given by 3 . 3-5. 

3.4 Star Tracking 

Unless a star tracker is a necessary part of the carrying 
vehicle's navigation equipment, it represents equipment provided 
solely for alignment. This implies that measurements of the 



dependent frame's inertial sensors are unable to provide a suitable 
alignment. The usual military application of star tracker align- 
ment relates to aligning a pre -leveled inertial system in azimuth; 
this requires one star sighting. Full three axis alignment is 
possible with two star sights. 

In spite of the numerous design options available in the 

25 

system application of star trackers , the basic problem remains 
one of accurately determining the line of sight (LOS) to a star. 

The vector input to the star tracking process is electromagnetic 
in nature, coincident with the position vector from the star to the 
system. The line of sight is concerned only with this vector's 
direction, It is possible, with an appropriate star and tracker, 
to determine an LOS within several arc seconds, even in daylight. 

The independent reference frame for star tracker alignment 
is the inertial frame. A star tracker can measure directly in 
this frame, thereby not requiring a master inertial system. With 
accurate time and position information, it is possible to translate 
the inertial frame measurements to an earth fixed reference frame. 

Figure 3-5 (See also Figure 2-27) casts the star tracking 
problem in the context of Chapter II. The angle a is measured 
from the star line-of-sight to the tracker-indicated LOS. The angle 
j3 accounts for those situations (See Section 3. 4. 1 below) where 
the transfer of star LOS information to the inertial sensor 
package produces an apparent deflection in the indicated line-of- 
sight. The total error, 7 , in measuring the star LOS is measured 
from the star LOS to the platform indicated LOS. 

3.4,1 Discussion of j 8 for Star Tracking 

Reference 25 informs us that military-type applications 
usually do not allow star trackers to be platform mounted. This 
is normally a mission oriented constraint required to provide the 
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Figure 3-5 Error Geometry Appropriate to Star Tracking 


tracker with an unobstructed view of the celestial sphere. As 
mentioned earlier, the technology is available (when the situation 
warrants) to provide values of a on the order of a few seconds of 
arc. Unless f3 is restrained to this same order of magnitude, it 
represents the limiting factor on star tracker azimuth. 

A simple way to visualize /3 is to consider the whole 
angle description of the rotation between the dependent frame and 
a coordinate frame fixed in the case of the tracker. The projection 
of this angle along the indicated line of sight is of no consequence; 
the perpendicular component is the angle fi. This presents us 
with a simple rule of thumb: when employing star trackers 
physically separated from the inertial sensors, the error due to 
this separation is reduced by placing the line of sight as nearly 
coincident with the principal flexure axis as possible. 
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ft may contain significant static, quasi-static and dynamic 
portions. The static portion is a result of inaccurate installation 
of the tracker and system. It exists when, say, an aircraft is 
on the ground or a ship is in port. Quasi-static contributions 
to j3 arise from approximately constant angular deflections of the 
carrying vehicle related to the average stress condition changes 
induced by motion of the vehicle (e.g. , aerodynamic loading of 
wings and pylons) and by environmental changes (e.g. , thermal 
gradient loading). The dynamic portion of ]3 is caused by excita- 
tion of the angular vibration modes of the craft. Total relative 
angular excursions between two points in the same vehicle may 
well be as large as 10° (Reference 3), 

Since /3 is strictly an angle transfer problem, it is 
theoretically possible to compensate the non-rigid transfer 
medium by means of optical comparison between the tracker and 
inertial system. This compensation is attractive but is not 
without its practical problems (See Section 3.8). 

In the absence of data which conclusively prove for a 
given application that j3 is either negligible or compensable by 
non-optical means, proposals to use a physically separate star 
tracker without optical compensation of flexural base motion 
should be viewed skeptically. When mission and economic con- 
straints permit, alignment accuracies will be enhanced by mount- 
ing the star tracker integral with the inertial sensors. It is 
wasteful to measure quantities better than you can transfer the 
information about those quantities 

3.4.2 Comment on 0 for Star Tracking 

The basic optics and detector of a star tracker are the 
fundamental physical limitations on determining aline- f-sight. 

26 

A frequently used figure of merit which describes this limitation 
is called the "tracking accuracy, " This accuracy number is 
identical with the parameter q-. 
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Rotations about a starline cannot be distinguished by the 
combination of optics and detector. On this basis, the conical, 
constant probability density locus interpretations used in deriving 
a (Section 2. 4. 2. 1) are a valid interpretation of fundamental 
star tracking limitations. The reader is cautioned, however, 
that the uniform statistics afforded by a vector viewpoint, when 
processed through angle transformations to obtain a component 
viewpoint, may lose the interpretation used in this thesis. 

3. 5 Fix Monitored Azimuth 

Fix monitored azimuth is a method of employing check- 
point data to align a pre-leveled platform in azimuth. In common 
with all usages of checkpoint information* (position reference) 
to reduce the effects of system and component error sources, 
fix monitored azimuth (FMA) is completely dependent on the 
accuracy of a model. The fundamental premise of FMA is that 
azimuth misalignment is not only the dominant source of position 
error during the time interval considered but is also nearly 
constant in this period. Given this model, the principle of FMA 
is exceedingly simple. 

A pre-leveled system, free from error in every respect 
except a fixed azimuth misalignment (0 & ), begins to navigate at 
station A (Figure 3-6) with the expectation of arriving at Station B. 

NOTE: 

g is normal 

into paper 
B 


Figure 3-6 Lateral Error Propagation for Fixed Azimuth 

Misalignment 

^Reference 5, pp 312 ff, presents a concise description of several 
uses of position references. 
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The range from A to B is R n . Due to 6 , lateral position error 

U cl 

propagates approximately as 


<5L = R0 

a 


3 . 5-1 


so that after navigating a distance Rq the lateral error magnitude 
at station B is R n e . By measuring 6L at B and knowing the range 
vector Rq, the azimuth error existing at A is calibrated and may 
be corrected. 

In the context of Chapter II, the lateral error vector 
is the input to the alignment process. From oL we infer the 
discrepancy between the nominal and actual navigation paths. The 
magnitude of 6L^ specifies the magnitude of 6 a (equivalent to the 
general error parameter a) while the direction of 6L^ relative to 
the range vector determines the sense of the required azimuth 
correction. The reference coordinate frame is that in which the 
checkpoints are known. For direct azimuth alignment with respect 
to the earth, the checkpoints are physical observables (e.g., terrain 
features, radio or radar beacons) of known location. There is 
no reason, however, why any point in space where the vehicle's 
position can be calibrated should not suffice for the FMA method. 

Measurements of <5L^ generally will be imperfect. This 
leaves a residual uncertainty in the corrected azimuth at each 
checkpoint. A powerful means of reducing certain measurement 
errors is to use multiple checkpoints. Let us now examine this. 


Allow errors in the measurement of <5L^ to be of the bias 
(a n ) and proportional (Aa. ) types. The residual azimuth error 
(A0 ). at the i checkpoint due to these error sources may be 

cl 1 

expressed as 


(A 0a )i - 


a o + Aa i 


6L i 


R. , . 
i-l, i 


3 . 5-2 
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Assuming that an azimuth misalignment estimate has been made 
at the previous (i-1) checkpoint, (A Q ). is actually the uncertainty 
in estimating the previous error (A© ). . , 6 l. may in fact, be 

cl l” X 1 

expressed as 


<5 L. 

l 


R. 


i-1, i 


(A e )• 
a l- 


3.5-3 


As long as the proportional measurement error is dominant 
( Aa^ > > a^) the azimuth uncertainty at the i checkpoint 

is written (from Equation 3. 5-2). 


Aa, 6L. 

<^ a >i r7 L 

ai K i-l,i 


3. 5-4 


Substituting 3. 5-3 into 3. 5-4 yields the simple recursion 
relationship 


< Ae a’i “ Aa l <A< Vi-l 3 - 5 ' 5 

Equation 3. 5-5 is very significant. It implies that every time 
we make a misalignment estimate, we obtain a fixed percentage 
reduction in the error in our previous best estimate (Aa^ < < 1). 
Iterating 3. 5-5 from the first checkpoint initial azimuth error 

— | +Wi 

(© a )i to the n checkpoint, the error in estimating ( 0 ^ 
depends only on (0 ), and the number of checkpoints (See Equation 
3. 5-6). 

< A 9 a>n= (^l )" 1 «Vl 3 ' 5 ' 6 

th 

The reader is reminded that in going from the first to the n 
point we have only estimated azimuth (n-1) times. 

The natural lower limit of this rapid convergence of the 
azimuth error is provided by the bias error (a Q ) existing in the 
measurement of <511^. When the lateral error is small enough 
so that bias becomes dominant, we see from Equation 3. 5-2 that the 
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lower limit on azimuth correction is expressed as 


(A Vi =~R 


i-1, i 


3.5-7 


Note that this is a function only of the bias measurement error 
and the range between the last two checkpoints. 

The convergence of the FMA process is so rapid that two 
or three inter-checkpoint legs suffice under most conditions to 
reach the bias limitation of 3. 5-7. In the face of a 10% propor- 
tional measurement error, a 10° initial azimuth error is* reduced 
to 6 arc-minutes after only two estimates. But on a 100 mile 
final leg, the bias error must be no greater than 150 feet to 
ensure azimuth alignment of 1 arc-minute. An arbitrary increase 
inR. . would seem to reduce (A 0 ). to an arbitrarily small 

1 “ 1 1 3 . 1 

level. The penalties here are two-fold. First, it is a severe 
operational restriction to require large inter- checkpoint distances. 
Second, as the time to complete the FMA process increases, the 
validity of the model deteriorates rapidly. 

Let us examine this latter point briefly. If a constant, 
uncalibrated platform drift (D o/hr) exists about the vertical, 
the lateral error due to this drift propagates as 

6L. s 1/2D t 2 3. 5-8 

1 v 


If a small constant platform drift (D R ) exists about the range 
vector R q , the lateral error propagates as 

6L. = 1/6 D R gt 3 3.5-9 

For demonstration purposes, we sum the contributions of 
3. 5-3, 3. 5-8, and 3. 5-9 to 6L^. 
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ir 

6L i “ R i-i, i < a v 1-1 + 1/2 v 2 + 1/6 “rS * 3 

3. 5-10 

If the vehicle maintains a roughly constant speed, V, the lateral 
error propagation due to azimuth error is proportional to the time 
elapsed between points (i-1) and (i). It should be clear that the 
model postulated for a system using FMA is valid for a restricted 
time interval; from 3. 5-10, azimuth error is dominant over the 
drift sources only if 

V(A 9 ). t. , . > 1/2 D t 2 + 1/6 g D t 3 . 
a 1- i 1- lj 1 V r\ 

The higher the speed of the vehicle, the more applicable the 

technique. Fix monitored azimuth is basically a fast vehicle 

technique suggested for airborne missile applications. 

Error sources other than system drift conspire to destroy 

the validity of the FMA model (e. g. , a steady state vertical error 

about R will cause an azimuth error to be introduced even if none 
o 

existed previously). The foregoing examples were merely by way 
of illustration. 

In this author's opinion the major deterrents to more wide- 
spread application of the fix monitored azimuth concept are 

1) The severe operational requirements imposed 
by navigating a required course 

2) The inherently long reaction time 

3) The restricted validity of the model on which 
FMA is based 

Mitigation of the first objection would be obtained by provision of 
a master inertial system whose continuous navigation outputs 
serve for checkpoint comparison at chosen intervals. This does 
not evade the other objections, however, leaving other master- 
slave configurations more desirable. 

We present here for the sake of completeness a variation 
on the FMA procedure. It is a proposed means for the azimuth 
alignment of a pre- leveled aircraft navigator/. Referring to Figure 3-6, 
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checkpoints A and B are at opposite ends of a calibrated runway. 

In other words, the direction of R is calibrated relative to true 

o 

north. If the pilot holds the course of the plane very close to R q 

during takeoff, then the average direction of the plane's horizontal 
acceleration is well-known. This is sensed by the navigator's ac- 
celerometers and used for alignment. 

There are three basic error sources in this scheme. 

First, if the pilot misses the takeoff point (B) by an amount 6L, 

the azimuth calibration is in error by 

0 = 6L/R . 

a o 

On a 1C, 000 foot runway this amounts to approximately 20 sec per 
foot of lateral error Secondly, even if the pilot passes directly 
over B, any lateral velocity is an indication that the average ac- 
celeration vector does not lie along R . Thirdly, there is the 

limited ability of the navigator's accelerometers to locate the 
average acceleration vector (given by the angle a of Chapter II). 

It is relatively easy to measure and compensate for the first two 
errors but we do not wish to consider this scheme in further 
detail here. As with the basic FMA technique, it suffers from 
limited applicability 

3. 6 Vector Matching 

3.6,1 " Memory" in Transfer Alignments 

Vector matching is the first transfer alignment 
technique (See Section 3, 1) we will consider. All transfer tech- 
niques are concerned with aligning a dependent reference frame 
to a pre-aligned intermediate frame. The inertial system asso- 
ciated with the dependent frame is often called a "slave", that 
associated with the intermediate frame a master". 

The analysis of transfer alignment errors relates 
to errors developed between the master and slave only. The mas- 
ter system is assumed perfectly pre-aligned to an independent 
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reference frame. Any error existing in the master's "memory" 
of the independent frame is necessarily transferred to the slave 
system. The total error in mechanizing a copy of the independent 
reference frame in a slave system is a linear superposition of the 
error in the master's memory and the error in aligning the slave 
to the master (See the functional diagram 2-1). A simple means 
of evaluating the total error for specific cases is by superposing 
the whole angle rotations representing master system drift and 
master-slave misalignment. 

3. 6. 2 Vector Matching Theory and Mechanizations 

The concept of vector matching is to constrain each 
of two systems to have an identical attitude relative to a common 
vector. Since the vector acts only as an intermediate variable, 
its orientation in the independent reference frame is not con- 
strained (and |3 is automatically zero). In the case of inertial 
vectors this means we can utilize "any" acceleration vector, not 
just g; we can use " any" angular velocity vector, not just u ^ . 

Gross motions of the carrying vehicle are turned to good advan- 
tage as generators of vector quantities measurable by the system's 
inertial sensors*; they are no longer interfering motions. Nomi- 
nally horizontal acceleration or velocity change vectors can be 
generated for azimuth alignment purposes by simple vehicle 
heading changes. 

Although the concept and basic mathematics of 
vector matching apply to every .measurable vector (inertial, 
magnetic, electromagnetic) we restrict out attention to those 
measurable by inertial sensors**. Earlier analyses distinguish 

27 

* Sutherland applies optimal linear estimation to maneuver 
determination for an airborne problem. 

** The special case of "vector matching" by means of a light beam 
is discussed in Section 3. 8. 
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between acceleration matching and velocity matching in a manner 
which implies the existence of fundamental differences. These 
differences are neither conceptual nor mathematical. The vec- 
tor V discussed in the following development represents any vec- 
tor, be it acceleration (linear or angular) or velocity (linear or 
angular). Any and all differences must lie in 

1) the ability to generate a particular vector 

2) the instruments available to measure 
vector quantities 

3) the resolution of the instruments. 

The mechanization of vector matching takes several 
forms. These are well illustrated by the two-dimensional problem 
of Figure 3-7, corresponding to the historically important case of two 
platforms which have already been erected to the vertical (i. e. , 

Z axes coincident) by some prior operation and which are to be 
accurately aligned in azimuth by means of a vector matching pro- 
cedure. If the two equivalent coordinate frames X^-Y^ and X 2 -Y 2 

are identically aligned to the vector V (but not necessarily about 
the vector*), the magnitudes of the components of the vector 
measured along corresponding instrument axes will be equal. 
Conversely, when the coordinate frames are misaligned (by the 
angle <j> x = <j> =$ ) the measured components will be different. 

Constraining the measured magnitude components to be equal by 
slewing the slave system results in proper alignment. This mech- 
anization seeks to zero the orthogonal magnitude differences: 

AV x = V xl - V x2 = V [cos ip - cos (0 + <l> )^j 

3. 6-la 

or 

AV y = V - V y2 = V [sin if, - sin (0 + cj> )] 

3.6-lb 

^Alignment about V is achieved by vertical erection. 
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Figure 3-7 Geometry for Azimuth Alignment 
by Vector Matching. 



t J 1 (R+AR) 
|w 2 R 


Figure 3-8 Different Magnitude, Parallel Direction 
Accelerations and Velocities Sensed by Two Systems 
on a Common Base 
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where V = [ V | . Applying Taylor f s series expansions about ip to 
Equation 3. 6-1 yields to first order in <f> 


or 


(j) 


AV /V AV 


x 


sin ip V 


yi 


3. 6-2a 


AV /V AV 
* 3 r -r =TT JL 

cos * V Xl 3. 6-2b 

Due to measurement errors, AV x and AV^. will generally have 

non-zero minima with a corresponding residual misalignment. 
To the extent that there are equal irreducible minima in the x 
and y channel measurements (i. e. AV^. = AV^ = AV), it becomes 

obvious that a judicious choice of the mechanized difference 
quantity guarantees a maximum azimuth misalignment of 


<f> = V2 

max V 


3. 6-3 

This occurs at ip = 45 u and agrees with a basic result of Chapter 
II. Equation 2.4-71 shows that independent of the vector match- 
ing mechanization, the error of each system in determining the 
direction of V is maximum at ip =45°. 


There may be situations where it is more advanta- 
geous to numerically evaluate the misalignment angle rather than 
eliminate it by a continuous difference technique. The computed 
angle may then either be (1) stored within the system computer 
and accounted for in all subsequent computations or (2) torqued 
out by a precision torquing procedure. Equation 3. 6-4 evaluates 
0 in terms of measured magnitude components. 





V 

-1 y 2 
-i p = tan -y — 

x 2 


V 


- tan 


-1 


3. 6-4 
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3. 6. 3 Compensation for Platform Separation 


It has been shown in Chapter II that alignment requires 
finding vector directions; magnitude measurements may be a 
means to that end (e. g. , Equations 3. 6-1 and 3. 6-4). A physical 
separation of two inertial systems attached to the same vehicle 
can cause both magnitude and direction differences in the maneuver- 
generated vector as sensed by each system. 

Consider as an example a vehicle traveling in a perfect 

circle (Figure 3-2). The magnitude of the centripetal acceleration 

sensed by two systems separated a radial distance All differs by 
2 

an amount ARw even though the sensed vectors are collinear. 
Similarly, the tangential velocity vectors differ in magnitude by 
ARu even though the sensed vectors are parallel. If the separa- 
tion were tangential rather than radial, and both systems were equi- 
distant from the center of rotation, then the magnitudes of the sensed 
vectors would be identical but their directions would differ. 

Directional differences must always be compensated; mag- 
nitude differences, when employing schemes such as 3. 6-1 and 
3. 6-4, must also be compensated. This is readily accomplished 
by measurements of the master system together with information 
on relative system locations. 

3. 6. 4 The Vector Match Vertical - A Fallacy 

The existence of a pre-aligned master system which can 
calibrate horizontal accelerations of a vehicle presents a strong 
temptation to establish vertical alignment by vector matching in 
the horizontal plane. The following statement should therefore 
be read carefully: It is impossible to achieve accurate vertical 
alignment by acceleration matching in the horizontal plane unless 
an accurate azimuth alignment has been established first. This 
is simply because a component of g projected on an axis due to 
vertical error is indistinguishable from a projected component of 
horizontal acceleration. Yet both are subject to the same 
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interpretation, namely vertical error. Only if azimuth alignment 
is perfect will the horizontal acceleration components sensed by 
the master and slave be equal and cancel when component dif- 
ferences are examined. 

The basic fallacy arises because we are dealing with com- 
ponents of a single vector. As discussed in Chapter II, sensors 
are only insensitive to rotations about the total vector. Rotations 
about one component are still sensitive to the orthogonal com- 
ponents, posing the problem of ambiguous data which must be re- 
solved by other measurements. When dealing with components of 
a vector, one must keep track of all the components at all times. 

Although vertical erection by vector matching in the null 
plane of g has been shown improper, it is correct to consider the 
normal vertical erection process as a vector matching process 
in that one attempts to match the best estimate of the position of 
the gravity vector to a pre- determined attitude. It differs from 
the ordinary concept of vector matching in two ways: 1) the vec- 

tor is naturally existing and 2) the desired relative orientation is 
known with respect to an independent reference frame, removing 
the need for a master system. 

It follows from this latter property that the possibility 
exists to require no data transfer link between master and slave 
platforms during a vector matching alignment, even though the 
vector is artificially induced. This requires that the orientation 
of the vector be tightly controlled to a pre- determined orienta- 
tion with respect to the master system reference axes. Since the 
vector is induced by a maneuver of the carrying vehicle, this is 
tantamount to precision control of the vehicle's dynamic path. A 
well known example of this type of control is maintaining the thrust 
vector of a missile through its center of gravity while adjusting 
the missile attitude with respect to the thrust vector. 


128 



The simplicity and reliability gained by use of this 
approach are due mainly to a reduction in data paths and system 
interconnections. These advantages come only at the high cost 
of an obsolute restriction on discretionary maneuvers of the 
vehicle pilot or controller during the alignment period. The 
accuracy of alignment obtained with communication between 
master and slave would be somewhat degraded by discretionary 
deviations from a prescribed maneuver pattern but this error 
would be bounded; without master-slave communication the error 
is obviously unbounded. 

3. 6. 5 General Formulation of Vector Matching 

A single vector match is capable of fixing two degrees of 
angular freedom. The solution of the azimuth alignment problem 
by vector matching, since it reduces only one degree of angular 
freedom, is a special case. We investigate here the more general 
vector matching problem. 

As shown previously, the inability of a system of sensors to 
perfectly indicate the direction of an existing physical vector can be 
resolved into a single normalized error angle, a. When the mea- 
surement apparatus is not directionally constrained about the 
measured vector, the constant probability density error locus is a 
cone of half-angle a(Figure 3-9a). The tip of the mis -indicated 
vector describes a cfircle in a plane normal to the measured vec- 
tor. All points on this circle have a constant probability density. 
Normalizing the measured vector sets the circle's radius at a. 

For separate systems using different instruments, the 
descriptions of the errors made in measuring the direction of 
the match vector have the identical geometrical interpretation 
and, insofar as the effects of environment can be ignored, are 
statistically independent. Assuming that the corrections have 
been made to ensure collinearity of the normalized match vector 
(Section 3. 6. 3) at all systems, the complete error analysis for 
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MIS-INDICATED 

DIRECTION 


(a) 


CONSTANT PROBABILITY DENSITY LOCUS: 



Figure 3-9 Error Geometry for Vector Matching 


n systems may be carried out in a plane diagram consisting of n 
concentric circles centered about the match vector. The case of 
greatest practical interest is n = 2, which describes the alignment 
of a single slave to the master. The resultant misalignment $ 
between the dependent and intermediate frames for this case is 
defined in Figure 3- 9b. The angle r? is identifiable as the difference 
of two independent angles, each uniformly distributed on 0-2n, 
Application of the law of cosines to Figure 3-9b yields: 


2 2 2 

= a + a - 2a a cos p 3. 6-5 

ms ms 

the statistics of which are evaluated in Appendix B with the results 


Mean Value: $ =— [ (or +a )£ (k 2 ) ] 

7 T m s 


3. 6-6 


2 2 2 

Mean Squared Value: $ = a +a 

ms 


3. 6-7 


where£(K ) is a complete elliptic integral of the second kind and 

3. 6-8 


k 2 =. 


4 a a 
m s 


( a + a ) 
m s 

Equation 3. 6-6 cannot be further evaluated without complete 
knowledge of the probability density functions for and a g . 

In keeping with the precedent established in Section 2. 4. 2. 1. , 
we will use the mean square misalignment as our performance 
index. 

The master system will ordinarily make more accurate 
measurements than the slave, i. e. if we make the identification 


2 2 2 q C Q 

a s q “m 3 * * - 6 " 9 

then q j^l. q may be thought of as a quality factor or figure of 

merit upon which systems are compared. It relates the RMS 

total vector measurement capability of systems of sensors as a 
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function of measurement method, number of sensors and quality 
of the sensors, q is very simply calculated by using Equation 2. 4-49, 
Tables 2-2 and 2-3, and a little algebra. The mean-square 
misalignment resulting from a vector match is expressed as a 
function of q and the measurement capability off the master system 
as 


$ 2 = (q 2 + 1) ; q >1 

or in terms of RMS quantities 


3. 6-10a 


RMS $ 
RMS a 


= (q 2 + 1) 1/2 ; q >1 


3. 6- 10b 


m 


Note, for example, that with identical instrument configurations 
in the master and slave but with an across the board order-of- 
magnitude saving in instrument quality (n = 10) in the slave, the 
contribution of the master system to the total misalignment is 
almost negligible. 


In order to consider full three degree-of-freedom align- 
ment by two vector matches, one merely obtains the angle <J» for 
each of the vector matches and uses the results of Section 2. 4. 3 
to calculate a composite misalignment for the entire system. 


3. 7 Gimbal Angle Matching 

The generic term "gimbal angle matching" implies the 
comparison of signals from the attitude readout systems of each 
of two (or more) inertial platforms in order to bring the orien- 
tations of their stable members into angular coincidence. Here, 
as in all methods involving comparison, the inertial orientation 
of one platform (master) is taken as an absolute reference and 
all errors in alignment of the second platform (slave) are with 
respect to this reference. Since the only controlled function is the 
orientation of each system's stable member with respect to the 
base on which it is mounted, there is a strong dependence on the 
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fidelity of reference information transmission by the base. This 
dependence is the most significant limitation of gimbal angle 
matching. Knowledge of this limitation is not new. The method 
of developing it here, however, is a new approch. 

Although a definite aid in the visualization of gimbal angle 
matching, there is no requirement that the systems have identical 
gimbal arrangements, nor must they even have an equal number of 
gimbals. Accurate computation of appropriate constant trigono- 
metric transformations of the attitude drive signals is technologi- 
cally feasible and can adequately compensate for this difference. 
The treatment of any desired non- coincident axis orientation 
follows similarly. An example of the latter situation is the align- 
ment of a preferred- axis -orientation missile platform to a carry- 
ing vehicle's geographically oriented navigational platform. 

Gimbal angle matching is the only alignment technique 
discussed in this thesis which is not applicable to "gimbal-less" 
or "strapdown" inertial measuring units. For what insights it 
may provide, however, the forthcoming analysis does apply to 
strapdown systems if the transformations T^ and Tg are replaced 
by the identity matrix. 

3.7.1 Coordinate Transformations 

Referring to Figure 3-10, an arbitrarily chosen set of 
orthogonal axes associated with the stable member of the de- 
pendent system is related to an orthogonal coordinate frame in 
its base by 



and further related to the coordinates at the base of the reference 
system by 
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MASTER 

REFERENCE SYSTEM DEPENDENT SYSTEM 



Figure 3-10 Symbolic Diagram of Gimbal Angle Matching Transformations 



r 


3. 7-2 

where Tg and Tg are square matrices of order three. For the 
master system we may write 


3. 7-3 


Equating 3. 7-2 and 3. 7-3 



. r l 





d l 

T 1 

r 2 

* 

T 

L 3 


T 2 

<N 


r 3_ 





d 3_ 


3. 7-4 

relates the stable member coordinates of the two systems. The 
elements of and Tg are functions of their respective system's 

attitude signals and possibly some constants; the elements of Tg 

are functions of the parameters chosen to describe flexure of the 
base. This flexure may be divided into quasi- static and dynamic 
parts. The quasi-static flexure is an approximately constant 
deformation of the base which varies with the average stress 
conditions induced by conditions of motion. The "droop" of an 
aircraft wing as a function of fuel load, speed, etc. is an example 
of quasi-static flexure. All other base flexure constitutes the 
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dynamic portion. In addition to parameters of flexure it is con- 
venient to include in the matrix Tg those constants necessary to 

describe the initial relative misalignment of the bases occurring 
at the time of system installation. 

Premultiplying both sides of 3.7-4 by the inverse of 

(denoted by *) yields the more useful and traditional repre- 
sentation of gimbal angle matching 


T 1 T 3 T 2 



3. 7-5 


where the unbracketed letters T. represent the previously 
introduced 3x3 matrices. Without loss of generality we may 
assume that the desired alignment corresponds to parallelism 
between (r r r^) and (d^, dg, dg) respectively and also 


tit tit 

between (r^, rg) and (dj, dg, dg). Any transformation re- 


quired to force the physical situation into agreement with this 
condition is constant, is assumed arbitrarily accurate, and 
hence does not enter the ensuing error analysis. 


3.7.2 Error Analysis 

The condition of angular coincidence between the two 
stable members is mathematically defined by 


Tj" 1 Tg T 2 = I 3. 7-6 

where I is the identity matrix. Under the simplifying restrictions 
of Section 3. 7. 1 and temporarily neglecting errors in the attitude 
systems we may specialize 3. 7-6 to require 
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3. 7-7a 


T 


1 


= T 


2 


T 3 “ 1 3. 7-7b 

With 3. 7-6 and 3. 7-7 as our definition of the desired alignment, 
we proceed to investigate the errors inherent in gimbal angle 
matching. The whole angle rotation concept of Chapter II will 
be useful here. 

It is convenient to define the two attitude systems by 

T 1 =T 1 + € 1 3. 7-8a 

/ 

T = T 4* e 

l 2 2 2 3 . 7 - 8b 

where a primed transformation corresponds to the ideal commanded 
attitude of a stable member and e accounts for the difference 
between the actual and commanded attitudes caused by imperfect 
attitude readout equipment and imperfect attitude drives. The 
"matched" condition 3. 7-7a is now T^ = T ^ and 3. 7-8b may be 
written as 


T 2 =T 1 + e 2 

or using 3. 7- 8a 


3. 7-9 


T 2 = T 1 +{e 2 6 1 ^ 
Substituting 3. 7-10 into 3. 7-5 gives 


3. 7-10 


T l~* T3 Tj + T i _1 T 3 (^-^>*^.2 


3. 7-11 


as the transformation between stable member coordinates when 
flexure of the base and an imperfect attitude system are considered. 
The transformation T^ ^ is constrained to be orthogonal. 


By Equation 2.4-15, the magnitude of the whole angle rota- 
tion described by is evaluated from the trace of Tj_g • 
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Since the trace of a sum of matrices equals the sum of the traces 
of each matrix, the error contribution of each term of 3. 7-11 to 
$ may be calculated separately. For the first term of 3. 7-11 

tr T 1 T 3 T l =trT 3 3.7-12 

This follows by application of the similarity argument developed 
in Equations 2.4-12 through 2.4-15. The physical implication of 
this result is that for assumed perfect attitude systems (e^ ~ e 2 = 0) 

the misalignment angle magnitude is only a function of the para- 
meters of the base and the initial installation. This rather simple 
(but important) insight is emphasized by using the whole angle 
rotation for the error parameter. Although it is frequently a 
matter of engineering practice to neglect attitude system errors 
in the face of much greater base flexure errors, the alignment 

3 

errors are still calculated as a function of the gimbal angles . 

This results from carrying out the multiplication T^ ^ T 2 to 

obtain the direction cosines relating stable member reference 
axes and a subsequent appropriate combination of these elements 
to obtain azimuth and vertical errors. All that is accomplished 
by this lengthy calculation is to find one particular coordinatized 
representation of the basic misalignment angle $ . Concern with 
the mathematical processes involved in multiplying three square 
third order matrices and the complexity of the resulting trigono- 
metric error expressions can easily obscure the simple geometry 
of the problem. The misalignment angle and the axis about which 
it is defined can be determined by reference to the base alone and 
then treated as geometric entities in space. The reduction of this 
error into components in a particular coordinate system then pro- 
ceeds directly without requiring reference to intermediate quantities. 
If the error axis is visualized as fixed in the carrying vehicle 
(base) and the attitude of the vehicle is specified with respect to 
a reference set of navigational coordinates then the decomposition 
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of the alignment error into components in the reference frame 
follows in a straightforward manner. This approach permits 
visualization of the change in alignment error as a function of 
vehicle attitude without resort to the substitution of "representa- 
tive" numbers into a complicated trigonometric expression in 
order to further understanding of the result. 

Perhaps more importantly, it suggests the possibility 
of dynamically compensating the errors due to installation and 
flexure. The measurement of inter-base angular phenomena 
could be achieved by an optical calibration link, such as a three- 
axis auto-collimation device which requires but a single light path. 
The cost-penalty of introducing equipment related only to align- 
ment is minimized by affixing the relatively inexpensive passive 
portion of such systems to the dependent system. The more 
expensive active mechanism remains with the carrying vehicle. 
Because the link is between bases fixed relative to the vehicle, 
turns of the vehicle do not disrupt the optical path. 


In the absence of a means for real-time base motion 
calibration, the accuracy realizable by gimbal angle matching 
is, for most vehicles, very crude (on the order of degrees). 

This technique then finds application only as a rapid, coarse 
alignment preceeding some other more accurate technique. The 
error bound for gimbal angle matching is then determined by the 
uncalibrated static and quasi- static portions of tr T^; the dynamic 
portion of tr T^ may be eliminated by filtering. 

If compensation of tr T^ is achieved to the order of a 
minute of arc or so, one must then undertake the alegebra asso- 
ciated with the attitude system contribution to the total misalign- 
ment, $ . From Equations 2. 4- 15. and 3. 7-11 the relationship is 


cos $ =1/2 


ti 


tr T 3 + tr T x 


-1 


T 3 * e 2 



3. 7-13 
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While evaluating tr T^ 1 (eg " involves considerably more 
algebra then calculating tr T^, it is only one third the work of 
calculating the triple matrix product 

3. 7. 3 An Observa tion Regarding System Installation 

An observation regarding system installation on board a 
vehicle is prompted by the whole angle treatment. First, there 
are certain guidance system designs which can tolerate a cruder 
alignment about one axis than about the others. Second, the axis 
of the whole angle misalignment will have a statistically preferred 
location relative to the flexural mode axes of the vehicle. If the 
system axis mentioned above and the predominant flexural mode 
axis are made parallel, the alignment problem is reduced. This 
may require, however, that the platform be torqued to a different 
orientation prior to inception of a mission. 

3. 8 Optical Alignment 

Optical alignment techniques are inherently very accurate, 
very rapid and simple in concept. They have proved their merit 
at fixed launch sites, particularly for azimuth alignment. In the 
moving base environment, they are an apparently attractive means 
for providing one, two or three axis alignment of a dependent frame 
to an intermediate frame (i. e. , slave to master). 

Optical comparison techniques circumvent the structural 
deflections of non-rigid bases. This is accomplished by establish- 
ing a light beam as a common reference between the systems. 
Provided a beam can be maintained between the master and slave 
systems, the systems need not even be carried aboard the same 
vehicle; the orientation of each is slaved to the electromagnetic 
vector. In truth, however, sufficient practical limitations exist 
that optical alignment on a moving base is normally dismissed as 
an impractical application. To this author's knowledge it has only 
been operationally employed in the submarine environment. 
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Not the least of the problem areas is the provision of a 
line-of-sight between the systems. Permanent obstacles require 
that the beam containing alignment information be directed around 
these obstacles. The use of prisms, mirrors or periscopes to 
direct the beam will reintroduce to some degree the very struc- 
tural deflections of the base we are seeking to circumvent. Ad- 
ditionally, even in the absence of structural deflection every 
reorientation of the beam is accomplished at the expense of 
angular error and light loss. Personnel and portable equipment 
must also be kept out of the beam's path. This might prove a 
very unwelcome restriction on an aircraft carrier's flight deck. 

As with all alignment schemes, operational requirements and 
restrictions are a major basis for the pre-selection of feasible 
techniques. 

The necessary provision of a suitable angular range is a 
significant restriction. A sufficient ’window" must exist to 
1) allow initial acquisition of the reflected beam and 2) maintain 
the beam, once acquired, in the face of angular vibrations of the 
vehicle. At least one application (Reference 3) forsees a 10° 
Master-Slave deflection. In an auto- collimating system, the 
angular range decreases as the separation of monitoring and 
monitored positions increases (See Figure 3-11). 

The feasibility of slaving one stable platform to another 

utilizing polarized light techniques has been demonstrated on a 
28 

limited scale . Althoughthis procedure replaces the imaging 

problem of auto-collimation with simple light gathering, it is 

29 

still restricted as to field of view. Both auto- collimating and 

2 8 

polarized systems now have the ability to transmit three axis 
information over a single light path. 

The heading of the vehicle carrying the two inertial systems 
must normally be restrained while employing optical alignments. 

As shown in Figure 3-12, unless the stable members of both the 
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Figure 3-11 Demonstrating the Angular Range Decrease with 
Linear Range Increase for a Fixed Autocollimator 
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master and slave are caged to the vehicle structure during turns, 
it is impossible to keep the two mirrors in view. The caged mode 
of operation is an extremely unlikely one for a master navigator. 

If we are considering strapdown systems, of course, no restric- 
tion on the heading of the vehicle is necessary. 

The accuracy achievable with optical slaving depends on 
the equipment selected, the maximum range, and the operating 
environment (windows, thermal gradients, vibration, etc. ). 
Accuracies generally range from a few seconds of arc to about 
one minute of arc per axis. Equivalent whole angle misalign- 
ments are calculable from the formulas of Section 2.4.1 when a 
comparison with other schemes is desired. 

3. 9 Summary 

The major accomplishments of this chapter are: 

1) Presentation of the basic theory of seven 
distinct alignment measurement techniques; 

2) Demonstration of the capacity of the unified 
alignment theory of Chapter II to properly 
include all the techniques; 

3) Exposition within the context of Chapter II 
of the fundamental sensor and base motion 
limitations on alignment accuracy. 

It has been shown for three of the direct measurement 
techniques (vertical indication, gyrocompassing, and star track- 
ing) that base motion is responsible for an apparent deflection of 
the measured vector relative to the independent frame. Those 
filtering and compensation techniques required to restore a 
meaningful direction to the measured vector have been enumerated. 

The fourth direct measurement technique (fix monitored 
azimuth) has been shown to successfully circumvent the problem 
of vector deflection by choosing to measure an earth referenced 
position vector. However, operational constraints, the time- 
limited validity of its basic model and its reliance on special 
equipment conspire to reduce its applicability. 
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The transfer techniques have been shown to circumvent 
the vector deflection problem by employing a master reference 
system. The vector imputs, rather than possessing an a priori 
significance in an independent frame, are calibrated by the mas- 
ter system. 

In the case of vector matching, the instrument limitations 
have been shown to be of the same form as tor vertical indication 
or gyrocompass ing. The error magnitudes, however, vary with 
the amplitude of the match vector; the selection of a technique 
has set /3 to zero. 

Angular deflections of the vehicle have been shown to 
limit the accuracy of gimbal angle matching. Although an old 
result, this was derived in a new fashion and provided additional 
insights regarding preferred installations. 

Optical comparison techniques have been shown to circum- 
vent the angular deflections of non-rigid vehicles. The only 
drawbacks to this technique, other than linear and angular range 
limitations, have been described as operational in nature. It 
remains, where feasible, an extremely attractive technique. 

There has been introduced in Section 3.6 an important, 
relative -quality factor, q. This factor describes the relative 
sensor limited capability of systems of sensors to indicate the 
direction of a common vector. That it can be easily calculated 
for sensors of varying performance, non-identical configurations 
and varying orientations of the sensor configurations relative 
to the vector input is a result of the consistent error format of 
Chapter II. 


144 



CHAPTER IV 


EXAMPLES OF UNIQUE OPERATIONAL SITUATIONS 


4. 1 Introduction 

Two examples of unique operational situations are dis- 
cussed in this chapter. The first case, that of simultaneously 
aligning several dependent systems in an airborne environment, 
reveals that measurements of a common input by multiple sys- 
tems promises no simple improvement in the alignment of these 
systems over that obtained independently. The second case, that 
of aligning aircraft navigation systems on a carrier at sea, demon- 
strates the degree to which operational considerations act as a 
pre-filter on any alignment procedure selection. 

4. 2 Multi ple Airborne Systems 

A useful concept in weapons systems technology is that of 
a highly maneuverable aircraft capable of employing multiple mis- 
siles. The independent alignment of each missile's guidance sys- 
tem relative to a master system or to an independent reference 
frame follows straightforwardly from the techniques of the pre- 
ceeding chapter. We now inquire whether the simultaneous align- 
ment of several systems is advantageous from the viewpoint of 
improving the alignment accuracy of each individual system. 

This question is most relevant to the simultaneous align- 
ment of n systems by measurement of a common physical vector, 

V (See Fig. 4-1). It is perfectly true that thip provides n indepen- 
dent measurements of a single quantity. It is also well known 
that multiple measurements of a scalar quantity reduce the effect 
of noise (or error) in the measurement process. To the extent 
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Figure 4- 1 


Measurement by n Systems of One Vector 


that the errors in each measurement are independent and that 
their variances are equal, (i.e. "identical systems") the va- 
riance in the sample mean of a scalar quantity is reduced by 

* /n. Let us see how this applies to the case at hand, namely 
multiple measurements of a vector quantity. 

The important parameters of V are its magnitude, V, 
and its direction. V is not only scalar, it is invariant under 
coordinate transformations. Therefore, the measurements of 

V by each system represent independent estimates of the same 
quantity. Simple averaging provides an improved estimate of 

V which is meaningful to each system; the magnitude estimation 
capability of each individual system, however, has not been im- 
proved at all. Overall knowledge of V is improved only because 
the invariance of V under coordinate transformation makes the 
measurement of V independent of the frame of reference. 

The direction of V must be estimated in order to determine 
alignment. Although the direction of V is representable by scalar 
quantities, these scalars (as opposed to V) are not invariant under 
a coordinate transformation; direction (i. e. , alignment) is a rela- 
tive quantity. Different systems' measurements of the direction 
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of V are statistically independent estimates of unrelated varia- 
bles. The fact that a common vector is being measured fails, 
in this case, to establish any dependency among the measure- 
ments. In the absence of additional measurements to establish 
some dependence among the variables, multiple alignment does 
not hold forth the promise of a simply achieved improvement in 
alignment accuracy (as by averaging measurements or "voting"). 

There is yet another reason why additional measurements 
are required to improve alignment accuracy. In the measure- 
ment of V, the quantity of interest (magnitude) is a physical ob- 
servable. In the measurement of the direction of V, however, 
the misalignment angle a is not a physical observable. Thus, 
the additional measurements required to improve alignment ac- 
curacy not only introduce some dependence among the variables 
but also serve the more basic purpose of rendering a physically 
observable. 

In answering whether the simultaneous alignment of se- 
veral systems is advantageous from the viewpoint of improved 
alignment accuracy for each system, we summarize here as 
follows : 

1) The existence of a single vector which is measured 
by several systems allows us to improve easily our knowledge 
of the vector's magnitude but not of its direction; hence, align- 
ment is not automatically improved by the existence of multiple 
systems measuring a common input. 

2) Misalignments can be reduced by means of additional 
measurements . 

The available magnitude improvement suggests the possi- 
bility of performing an on-board calibration of certain sensor co- 
efficients. (There is a brief calibration discussion in Appendix 
A. ) An improved knowledge of sensor coefficients, of course, im- 
plies an improved alignment capability. The precise feasibility 
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of this calibration and the operational restrictions required to 
achieve it will be revealed only by a detailed analysis of speci- 
fic situations. This is not considered further here. 

Let us consider, by means of a simplified example, an 
additional measurement which does make it possible to decrease 
the alignment uncertainty of each system. (Because each missile 
is likely to be separately targeted, an improved composite esti- 
mate of misalignment is valueless unless it provides the means to 
upgrade each system's alignment. ) Consider the plane case of 
Figure 4^2, where each system's misalignment is described by a ^. 

Although is not a physical observable, we wish to reduce 

its variance by measurements additional to those of the vector. 

The additional measurement we choose is the 9. , . between two 

systems' independent indications of V. Note that if n is even, 
only n/2 additional measurements are implied. These measure- 
ments are between missile systems, not between the missile and 
a master or missile and the independent frame. When the mis- 
siles are close together, as in a magazine, this angle may be re- 
latively easy to measure. We consider only the pair of systems 
shown in Figure 4-3 and assume the and have a known cova- 
riance matrix and a zero mean. The measurement of 9 is im- 
perfect, being corrupted by additive noise v of zero mean value 

2 

and known variance c . v is independent of and . 


Prior to measuring 9, our best estimates of and 

are zero (their means) with errors in these estimates whose 
2 2 

variances are a and a (the variances of the alignment er- 
a l a 2 

rors). After we measure 9 and create new estimates a ^ and a 2 
of aj and eg , the variances of the errors in our estimates of the 
misalignments are given by the results of Appendix D as 
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Figure 4-2 Definition of the Measured Angles 0^ i 



Figure 4-3 Definitions of Angles for Example of One System- Pair 




The measurement of 9 is seen to provide an improved estimate 

of the misalignments of each system relative to the vector. 

The actual estimator for is found in Reference 4. We have 

merely shown here the reduced variances of the errors in the 

estimates of misalignment provided by a measurement of 9. 

2 

Note that if the measurement noise variance a becomes too 

large, the measurement loses its effect. We plot D-9 and D-10 

2 2 2 

in Figure 4-4 for the special case a = a = a . For systems 

1 2 

with equal variances the largest possible improvement in the va- 
riance of t* is seen to be a factor of two. Whether or not this im- 
provement warrants implementation of the measurement of 9 and 
the computational routine which derives ^ and is a judgment 
which cannot be made here. 

In summary, this section has shown that multiple inde- 
pendent measurements of a single vector allow an improved es- 
timate of the magnitude of the vector to be simply obtained. 

These measurements do not provide a corresponding simple im- 
provement in the alignment of each measuring system relative 
to the vector. Additional measurements which relate the systems 
by pairs do make possible a significant alignment improvement 
for each system relative to the vector. 

4. 3 A Naval Situation 

The Navy's constant emphasis on reducing pre-flight 
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alignment times for carrier-based aircraft has resulted in two 
rather operationally unique alignment techniques. These tech- 
niques are known as the "Transfer Align" and the "Insertion 
Method. ' 

The present alignment of an aircraft inertial system at 
sea is a two-phase procedure; first, the system must be warmed 
to its operating temperature; and secondly, the alignment is 
achieved by gyrocompassing (which here means tracking the ho- 
rizontal component of earth rate from a pre- leveled platform; 
azimuth coupling is achieved via small platform tilts). The to- 
tal time required for the two-step process is the "ready time. " 

Personnel of the U. S, Naval Avionics Facility (Indianapolis) 

30 

have expressed the opionion that the ready time for gyrocom- 
passing in a parked aircraft "probably will never be reduced to 
less than 15 or 20 minutes. " At least half of this time period is 
devoted to gyrocompassing, during which the aircraft must be 
parked on the deck with uninterrupted power and reference signals. 
The length of time that the aircraft cannot be moved is considered 
a serious deck problem. 

The insertion technique was proposed as a means of re- 
ducing the overall ready time. In this scheme, a thermally sta- 
bilized, pre-aligned system is brought from some remote loca- 

30 31 

tion and inserted in the aircraft J ; the dependent frame, in 
other words, is portable. Rather than aligning the dependent 
system on a time- constrained basis, it has been continuously 
gyrocompassing or navigating and its alignment has been carefully 
trimmed by comparison with SINS. (SINS is the intermediate 
reference frame aboard naval vessels. ) The problem of aligning 
in a short time in the face of disturbing base motions has been 
traded for the strictly operational problems of hand carrying a 
system to the flight line and installing it in an aircraft without 
introducing severe thermal, mechanical or electrical transients. 
That this latter procedure can be consistently accomplished in a 
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short time has not, to this author's knowledge, been proved. It 
is, nonetheless, a very attractive possibility. Removal of the 
time constraint on alignment should permit greater gyrocompass- 
ing alignment accuracies with the same basic equipment. 

The transfer align method reduces only the alignment time; 

the permanently installed aircraft system must still be brought to 

operating temperature. This method employs an active transfer 

reference frame (Section 2.2.4) which is pre- aligned relative to 

SINS by-the techniques of Chapter III. The contractual transfer 
32 

align effort in support of the F-111B aircraft inertial naviga- 
tion system has chosen to use self-gyrocompassing of the trans- 
fer system with the assistance of SINS data. The transfer system 
is, therefore, an intermediary between SINS and the aircraft sys- 
tem; it may be thought of as an additional master system. Once 
aligned, the transfer system is carried to the aircraft (operating 
as a precision navigator on its own battery power) and attached 
to an external reference surface. Once mechanically attached, 

the alignment of the aircraft system is accomplished in less than 
30 32 

one minute J by means of a technique such as gimbal angle 
matching. 

In both the insertion and transfer align systems, the basic 
means of transferring alignment from one system to another 
(See Figure 2-1) is still chosen from those enumerated in Chapter 
III. But the emphasis is heavily on time. Time and accuracy 
are closely intertwined in this gyrocompassing mechanization, 
whereas the time scale is greatly reduced with several other 
alignment techniques. 

The insertion and transfer techniques attack the time pro- 
blem in distinctly different ways. The insertion method makes 
the time available for slow gyrocompassing by aligning completely 
"off line. " The transfer technique uses an inherently rapid align- 
ment transfer "on line. " The rapid techniques (i. e. , gimbal angle 
matching and optical comparison) require proximity between the 
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master and slave to eliminate base flexure and/or range problems. 
A preference for one technique over another is more likely to rest 
on operational considerations than on purely technical merit. 
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CHARTER V 


SUMMARY AND RECOMMENDATIONS 


5 . 1 Summary 

This study has attempted to develop a unified analysis 
of the alignment problem, with particular emphasis on the 
problems related to alignment on a moving base. The analysis 
facilitates a comparison and evaluation of the fundamental 
limitations of seemingly distinct techniques. It has described 
all alignment measurements as measurements of a vector 
quantity. Furthermore, it has demonstrated significant 
distinctions among "measuring a vector, " "measuring a vector's 
magnitude, " and "measuring a vector's direction, " where only 
the latter is of ultimate significance to the alignment problem. 
The fundamental inability of a system of sensors to indicate the 
direction of a vector has been reduced in all cases to an equiv- 
alent angular uncertainty. This uncertainty was shown to be a 
useful figure- of- merit for comparing systems of sensors. 

Angular deviations of the measured vector for a desired nominal 
direction have been ascribed to vehicle motions relative to the 
earth and flexibility of the vehicle structure. The combined 
effects of sensor imperfection and base motion represent the 
fundamental limitations on moving base alignment. 

It was desired that the basic figure- of- merit chosen 
for evaluating misalignments be insensitive to coordinate frame 



transformation, thereby rendering system-level comparisons 
independent of system geometry. The magnitude of the whole- 
angle rotation between coordinate frames met this requirement 
and also enjoys other important properties. For example, the 
magnitude of a rotation specified by the arguments of the princi- 
pal direction cosines is independent of the skewness of the 
coordinate frame, and the magnitude of a rotation specified by 
Euler angles is almost completely independent of the order 
of the rotations. The whole- angle concept made feasible a 
simply expressed evaluation of the sensitivity of misalignment 
to the angle between the two required vector inputs. In order to 
facilitate the adaptation of results expressed in terms of more 
familiar rotational parameters, the whole angle magnitude was 
developed in terms of many of these parameters. 

The importance of distinguishing between magnitude 
measurements and directional measurements has been demon- 
strated in the case of multiple system alignment. Simultaneous 
measurements of a single vector by multiple systems were shown 
not to constitute multiple measurements of a common alignment 
parameter. Therefore, no simple improvement in alignment 
accrued from aligning the systems simultaneously vs. aligning 
them separately. It was further shown that any misalignment 
resulting from only one vector measurement is not physically 
observable. Viewing the alignment problem as a vector-direction- 
indication process was crucial to these deductions and to 
suggesting an additional measurement which both renders the 
misalignment observable and introduces a useful redundancy. 

The relative advantage of this measurement was calculated. 
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The utility of the unified theory of alignment as an 
analytical tool was demonstrated by application to seven distinct 
alignment measurement techniques drawn from the literature. 

In addition to providing a common basis for comparing these 
techniques, it has been shown to be a simpler way of treating 
old problems (particularly gimbal angle matching). It also 
serves an important purpose beyond analysis. It enumerates in 
detail the basic requirements for alignment, and does this in a 
standard form. A simple test of sufficiency is thereby available 
prior to any detailed analysis of a newly discovered technique. 

5. 2 Recommendations 

The non-dimensional analysis of vector measurements 

has resulted in definite guidelines to minimize magnitude- indication 

and direction- indication errors as a function of instrument 

uncertainties and the orientation of the vector relative to the 

instruments. The results obtained are immediately applicable, 

yet this appears to be a profitable area for further study. 

Consideration of higher order instrument uncertainties is 

certainly appropriate, as well as whether redundant instruments 

(i. e. , a fourth gyro or accelerometer) would contribute 

significantly to a system's alignment capability. Studies such 

2 0 

as this have been made where the performance index 
relates to terminal navigation error, but not, so far as this 
author knows, where misalignment is the figure- of- merit. 

Another area of possible study relates to the 
calibration problem. In a crude sense, at least, the calibration 
problem may be viewed as the "inverse" of the alignment problem. 
That is, for alignment it is assumed that one has the best possible 
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knowledge of his instruments and seeks to ascertain the direction 
of a vector (sometimes with the help of magnitude measurements). 
For calibration, on the other hand, it is assumed that one has 
excellent knowledge of a vector and seeks to ascertain by means 
of this vector a detailed knowledge of his instruments. For 
those cases where the misalignment angle is a physical observ- 
able during calibration, certain geometrical and analytical formu- 
lations of this thesis carry over directly to the calibration problem. 
The extent to which these formulations aid an understanding of 
the calibration problem is of interest. Another calibration 
study was suggested in Chapter IV. It was shown there that it 
is relatively simple to obtain improved knowledge of a vector's 
magnitude when the vector is measured simultaneously by several 
systems. The possibility of using this information to update one's 
knowledge of certain sensor coefficients appears to this author to 
be worthy of investigation. The degree to which this updating 
improves one's ability to align each system is the performance 
index of interest. 

Finally, the author adds his voice to those who have 
previously emphasized the need for the improved availability of 
data concerning vehicle motions and descriptive models of this 
motion. The first requirement for a detailed dynamic study of 
any of the techniques described in this thesis is a set of 
representative models for motion of the various vehicles relative 
to the earth and for the relative linear and angular structural 
vibrations between inertial system locations. 
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APPENDIX A 


Instrument Example 
The Linear Accelerometer 


To clarify the definition of various instrument uncertain- 
ties, consider a simple linear accelerometer. 
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Point X is the center of mass of the "proof mass. " Point P is 
fixed in the case. The proof mass is mounted so as to have only 
one degree of translational freedom. 1 is a unit vector parallel 
to the axis of translational freedom and defines the sensitive axis 
of the accelerometer. The output indication of the instrument, x^, 
is made proportional to the displacement of the proof mass rela- 
tive to the case along 1 

S 3. 

Let m = the mass of the proof mass. 

k = the elastic coefficient of the spring, 
c = the damping coefficient. 

f = the non-field (i. e. contact) specific force acting on the 
instrument. 



p = the derivative operator 


d_ 

dt 


and write the governing equation 

z 2 . c , k x T • T 

(p + — p+ — ) x. = -f 1 
^ m ^ m 1 sa 


In the steady state for constant f 

x. = ~ ( - f * f ) + x. (0) 
l k sa i 

where x.. (0) is the zero setting of the scale. Choosing the 
nominal value x.. (0) = 0 permits x.. (0) to represent directly any 
uncertainty in the scale zero. If the elastic coefficient of the 
spring should be capable of randomly varying by an amount Ak 
but the mass m remains well known. 

x. = - ^ ( 1 - ~ ) f ' “ + x. (0) 

i k k sa i 


The foregoing equation describes an imperfect sensing mecha- 
nism which converts an input quantity of one set of dimensions 
(force per unit mass) to another (displacement). The uncertain- 
ties in this conversion are inherent to the transducer action 
referred to l ga - Hence the nomenclature "transducer uncertainty" 
to denote imperfections in knowledge of the basic conversion 
process. 


Normalizing the steady state equation for unity gain 

x.* = — x. = + — x. (0) 

l mi k sa m i 

and comparing with Equation 2. 4-55, the following identifications 

may be made 


0 


bias = — x. (0) 
m i 


Ak 


a. = scale factor = 1 - , 

1 k 


Ak 


Aa^= scale factor uncertainty = — ^ 


Under automatic operation which, for example, constrains 
1 flin order to place 1 parallel to f, the transducer 


sa 



uncertainties a„ and A a, will cause the angle between f and 1 

0 1 & sa 

to be non-zero. 

Sensitive axis uncertainty problems arise simply because 
the accelerometer is typically a sealed unit. This requires 
knowledge of 1 to be transferred to the case exterior as, for 
example, a scribed line. If the scribed line and 1 are skew. 

Set 

the sensitive axis uncertainty numerically equals the skew angle. 

Under automatic operation with no transducer uncertainty, 

1 can be made exactly parallel with the vector input f. To an 
external observer, however, the indicated vector direction is 
given by the direction of the scribed line. When the external 
observer role is assumed by the platform upon which the ac- 
celerometer is mounted, the platform becomes aligned imper- 
fectly in spite of perfect transducer behavior. 

In summary, sensitive axis uncertainty is measured from 
1 to the scribe mark while transducer uncertainty is measured 
from the vector input, f, to 1 . The total indicated misalignment 

S3 

is measured from f to the scribe mark. 


1 sa 




'sa 


SCRIBE MARK 


at 


a = total angular indication error. 

a = angular indication error contribution of sensitive 

S3 

axis uncertainty. 

a ^ = angular indication error contribution of transducer 

uncertainty. 
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To further distinguish between transducer uncertainty and 
sensitive axis uncertainty, consider the calibration of the linear 
accelerometer. The determination of bias and scale factor re- 
quires precise knowledge of the magnitude of g but relatively crude 
knowledge of its direction. Conversely, determination of the lo- 
cation of 1 requires precise knowledge of the direction of g but 

Set r 

relatively crude knowledge of its magnitude. To fully calibrate 
an instrument, then, requires complete knowledge of the magni- 
tude and direction of the vector used for calibration. Conversely, 
aligning an instrument with respect to a vector requires a complete 
knowledge of the magnitude and direction characteristics of the 
instrument . This contrast is the basis for viewing the calibra- 
tion problem as the inverse of the alignment problem. 


The determination of ag and a^ involves placing 1 both 
parallel and anti-parallel to "g" and taking respective readings 
(x.*) 1 and (x.*) 0 . Because the magnitude of g sensed by the ac- 
celerometer varies as the cosine of the angle between g and 1 , 

an uncertainty in this angle of up to 1 milliradian still guarantees 
knowledge of the input magnitude to at least five parts in 10 . a^ 
is found from 


a 


0 


2 


and a^ from 


a l 




- (x^*) 


2 



Thus, very accurate transducer calibration is had in the face of 
crude angular alignment. 

Having a n , the precise direction of 1 is found by rotating 

1 until it is normal to g (minimum indicated output). For an 
sa -7 I — I 

uncertainty (1 a) of 5 x 10 |g jin x * a Q has an uncertainty 
( 1 a) of 3. 5 x 10 ^ g. 
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This implies that 1 may be placed normal to g within 
— ^ sa 

3. 5 x 10 milli radian (or somewhat less than 0. 1 arc sec). 

An orthogonality accuracy of 1 milli radian can be had for a 

3 _ 

magnitude uncertainty of only 1 part in 10 . Placing 1 normal 
to g" within such tolerances is useful only if the direction of g (or 
its normal plane) is known at least as accurately, for only by 
reference to this direction can 1 be located meaningfully. 
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APPENDIX B 


A CERTAIN EVALUATION OF THE EQUATION 
A 2 = B 2 + C 2 ± 2 BC cos r? 


From the equation 

A 2 = B 2 + C 2 i 2 BC cos r] (B-l) 

under the conditions that r) is statistically independent of B and C 
and further that 

f] = ri 1 - r) 2 (B-2) 


where and rjg are independent variables, each uniformly dis- 
tributed on the int erva l 0 - 27T, we wish to evaluate the mean (A) 

2 

and mean- square (A ) values of A for B and C fixed. 


The probability density function for ri is given by 


f(r?) = 


47T 


~ (2 n + ri) 


Air 


^2 (2 7T - T]) 


-27 r < n <0 


0 < r) < 2 7r 


(B-3) 


The mean of A is derived from 
- ?' 

A ~ j A f (r]) d r] 

-27 r 

27T 


■ U 

- 2tt V 


>2 , ~,2 


1 


B + C ± 2 BC cos r) f (rj) d r\ 


(B-4) 


165 



In order to obtain an integrable form of (B-4), the integral is first 
treated in the following manner 

_ p 2 7T 1 

A = \ \ 2tt /B 2 + C 2 ± 2 BC cos n d 77 

l-i J 


o 




B 2 + C 2 ±2BC cos 


2tt 

h dr} - C r]I 
0 ' 


B 2 + C 2 ± 2 BC cos r) dn>(B- 5) 


Suitable changes of variable reduce (B-5) to (B-6). 

2n 2 7T 

A = j J ^B 2 + C 2 ± 2 BC cos nd r) - J j/b 2 + C 2 ± 2 BC cos rj 


dn 

(B-6) 


If (B-4) is treated as 

0 

1 J 
A 


• * { / 


B 2 + C 2 ± 2 BC cos rj (2 tt + n ) d r) 


27 r 


+ \ /b 2 + c 2 i 2 BC cos r] (2ir - rj) d rj 


(B-7) 


0 


the substitutions (27 r+r/) = 77 ’ and (2? r- 77 ) = 77 " in the first and second 
integrals of (B-7) respectively yield 
27T 


A = 


— 2 f rl/l 

2, 2 J / 


B 2 + C 2 ± 2 BC cos 77 d n 


(B-8) 


Equating (B-6) and (B-8) requires 
27T 27T 


^B 2 + C 2 ± 2 BC cos rjd n = jj 77 /b 2 + C 2 ± 2 BC cos n'dn (B- 

0 o'' 


9) 


which by comparison with the right hand side of (B-8) gives the 
desired result 
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(B- 10) 


2tt 





B 2 + C 2 ± 2 BC cos r) d n 


The form of the integral (B-10) is elliptic. By appropriate 
changes of variable (B-10) can be reduced to a tabulated form. 


For an example, consider the case where the plus sign 
occurs under the radical. The integral may be split as follows 

I ' * 

A = J 


2*' 


u 


1 


B + C +2 BC cos rj d rj 


(B-ll ) 


2tt 
7 r v 


B 2 + C 2 + 2BC cos T) d n 


Operate on the second integral of (B-ll) by first making the change 
of variable rj = (2n - n') and then interchanging the limits of 
integration. The two integrals of (B-ll) are then identical with 
the result 

~2 T— 1 

B^ + C + 2 BC cos n d r) (B-12) 

Equation (B-12) also results from manipulations on (B-10) when 
the minus sign occurs under the radical. 

If the substitution /3 = rj/2 is made in (B-12), the tabulated 

form of (B-13) is obtained. 

Tt / 2 

_ 9 p f 9 9 I 

A = — (B+C) j/l-k sin j3 d|3 

ov 



= V (B + C) £ (k 2 ) 
7 r 


(B- 13) 
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1 7 

£ is the complete elliptic integral of the second kind with 
argument k 2 


2 _ 4 BC 

(B + C ) 2 


(B-14) 


The evaluation of the mean- square value of A follows much 
more directly. A is defined by 

2ir 

2 f 2 / 

A = \ A f ( 17 ) d r) 

-2ir 


2n 

= j (B 2 + C 2 ± 2 BC cos 7 ])f (n) d r 7 (B-15) 

-2jt 


which is simply integrated with the result 

A 2 * B 2 + C 2 (B-16) 

The results (B-ll) and (B-14) are for B and C held constant. 

— 2 

If they are random variables, the dependence of A and A on the 
statistics of B and C would be calculated from (B-ll) or (B-14) 
respectively. 
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APPENDIX C 


AN INTERPRETATION OF STATISTICS ASSOCIATED WITH 
A TRIAD OF SENSORS 


Relative to measuring a vector with a triad of similar 
instruments, we may make the interpretation of Figure C-l: 
instrument uncertainties result in a three-dimensional uncer- 
tainty in the location of the tip of the vector V. The statistics of 
this uncertainty relative to the tip of the vector are obtained by 
affixing the x, y, z coordinate set to the tip of the vector such 
that it is a translated (but not rotated) representation of the x', 
y 1 , z' frame. If we now assume for convenience that the instru- 
ment uncertainties of Equation 2. 4-60 are normally distributed 
with zero mean, the statistics along x, y, and z are normally 
distributed with zero mean and standard deviations 


% ' J a B. 2 0 * y2 cos2 +1 


r 2 2 ~ 1 

a = / a 2 + V cos tli „ cr A 2 

y J a 0 2 Aa x 


CT z = 1 ° ag + y2 cos2 + 3 a A a 2 


(C- 


Note that although the statistics of each instrument are identical, 
the statistics of each axis are generally non- identical. This 
phenomenon results from scaling of the scale factor uncertainty 
contribution by the projected magnitude of the measured vector. 




I 

It is of great interest to many people to plot contours of 
constant probability density in order to exploit certain geometri- 
cal interpretations of probability. Unfortunately, this problem 
does not lend itself to this technique. The contours of constant 
probability density in three-space are described by 





Const. 


(C-2) 


The solid body described by C-2 is an ellipsoid, but one whose 
shape varies as its center is moved in three space. It becomes 
an ellipsoid of revolution if and only if two of the three direction 
cosines of the vector are equal. It becomes spherical if and 
only if all three direction cosines are equal. This, of course, 
only occurs when the vector is geometrically centered in the 
octant. There are exactly eight positions of the vector in all of 
three- space for which the contour of constant probability density 
is spherical! It should be obvious that for the majority of possible 
vector orientations with respect to the triad, the contours are not 
simple surfaces of revolution. 


Results such as 2. 4-66 are obtained without full calculation 
of the three dimensional error probability density functions for 
all vector orientations. The approach taken there is simpler and 
more direct. If one were to develop the three-dimensional pro- 
bability density and evaluate the statistics of a vector tip's loca- 
tion, great care would have to be exercised to see that the result 
obtained related to directional uncertainties and not just the 
"total uncertainty" in measuring a vector. 
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APPENDIX D 


OPTIMAL ESTIMATION APPLIED TO MULTIPLE SYSTEMS 


Consider two independent random variables a , and a 0 , 

2 2 “ 

each of zero mean and respective variances a and cr 




By one measurement involving a ^ and a ^ we can make better 
estimates'^ and of a ^ and a ^ such that our estimation error 
is improved 


(a ^ - &]) 




(D-l) 


2./\ . 

ct (a 2 ) 


< 


2 

<7 


a 2 


even in the face of measurement noise. This is done here by a 

straightforward application of optimal estimation theory. We 

9 

use here the formulation and notation of Bryson . 

The physical measurement is 

z = 0 + v = a 2 - + v (D-2) 


= H x + v (D-3) 

where the noise in the measurement is v and we have defined 


H = 


Pi, +ll , X = 

a i 

L_ 9 -J 9 

_ a 2 _ 


(D-4) 


v is assumed independent of x. with zero mean and known variance 
o 

c . Using the known relation of 6 to x. (Equation D-3), the known 
statistics of v, an initial estimate of_x, and the known uncertainty 
in this initial estimate, one can obtain an improved (in the least 
squares sense) estimate of _x (Reference 4). 
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The covariance matrix of the error in the estimate of x 
is given by 

P = E (e, e^) = E j(x - x) (x - x)^J (D-5) 

where Eis the expected value operation. P is given by 

P = M-MH T (H M H T + R)" 1 HM (D-6) 

where M is the covariance matrix of our initial estimation error 

2 


M 


= E (x - x) (x - x)^] = 


a-i 


0 

2 

Cf 

a 


(D-7) 


and R is the covariance matrix of the measurement error. 
R = E (v 2 ) = c t 2 


(D-8) 


Straightforward algebra yields for the variances in the errors in 
the estimates of a-, and a c 


CT ('a^ - a l )~ 


a 



a 


1 + 


a 


O', 


(D- 9) 


and 


^2 " 01 2 * 


= a 


a 


1 + 


2 

7 

a r 


a tc 


(D-10) 


The results D-9 and D-10 fit the restrictions D-l. 
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